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Abstract. We generalize the analysis of [12] and develop a singular pseudodiffcrcntial calculus. 
The symbols that we consider do not satisfy the standard decay with respect to the frequency 
variables. We thus adopt a strategy based on the Calderon-Vaillancourt Theorem. The remainders 
in the symbolic calculus are bounded operators on L 2 , whose norm is measured with respect to some 
small parameter. Our main improvement with respect to [12] consists in showing a regularization 
effect for the remainders. Due to a nonstandard decay in the frequency variables, the regularization 
takes place in a scale of anisotropic, and singular, Sobolev spaces. Our analysis allows to extend the 
results of [12] on the existence of highly oscillatory solutions to nonlinear hyperbolic problems by 
dropping the compact support condition on the data. The results are also used in our companion 
work [6] to justify nonlinear geometric optics with boundary amplification, which corresponds to 
a more singular regime than the one considered in [12]. The analysis is carried out with either an 
additional real or periodic variable in order to cover problems for pulses or wavetrains in geometric 
optics. 
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1. Introduction 



Nonlinear geometric optics is devoted to the construction and asymptotic analysis of highly oscil- 
latory solutions to some partial differential equations. In the context of nonlinear hyperbolic partial 
differential equations, one of the main issues is to prove existence of a solution to the highly os- 
cillatory problem on a time interval that is independent of the (small) wavelength. Such uniform 
existence results cannot follow from a naive application of a standard existence result in some func- 
tional space, say a Sobolev space H a , because the sequence of data does not remain in a fixed ball 
of H s . A now classical procedure for proving uniform existence results is to work on singular prob- 
lems with additional variables and to prove uniform energy estimates with respect to the singular 
parameter. This strategy was used in [10] for the hyperbolic Cauchy problem and adapted in [12] 
to hyperbolic initial boundary value problems. While energy estimates in [10] relied on symmetry 
assumptions and integration by parts, those in [12] are much more difficult to obtain and rely on 
a suitable singular pscudodiffercntial calculus. The operators are pscudodifferential in the singular 
derivative d x + fldg/e. This calculus is well-adapted to boundary value problems that satisfy a 
maximal energy estimate, that is an L 2 estimate with no loss derivative. In particular, remainders 
in the calculus of [12] are bounded operators on L 2 whose norm is controlled with respect to some 
parameter 7. This calculus is adapted to the situation studied in [12] because such terms of order 
can be absorbed in the energy estimates by choosing 7 large enough. 

In [7], two of the authors have studied and justified geometric optics expansions with an ampli- 
fication phenomenon for a certain class of linear hyperbolic boundary value problems. For linear 
problems, uniform existence is no source for concern. In the companion article [(>], we extend the 
result of [7] to semilinear problems. One major issue in [6] is to prove that the amplification phe- 
nomenon exhibited in [7] combined with the nonlinearity of the zero order term does not rule out 
existence of a solution on a fixed time interval. Our strategy in [6] is to study a singular problem for 
which we need to prove uniform estimates. As in [7], the linearized problems in [G] satisfy a weak 
energy estimate with a loss of one tangential derivative. 1 Such estimates with a loss of derivative 
were originally proved in [5] and are optimal, as proved in [7]. The amplification of oscillations is 
more or less equivalent to the loss of derivatives in the estimates. Compared with [12], we now need 
to control our remainders by showing that they are smoothing operators, otherwise there will be no 
way to absorb these errors in the energy estimates. Moreover, since the nonlinear problems of [G] 
are solved by a Nash-Moser procedure where we use smoothing operators, it is crucial to extend all 
the results on the singular calculus of [12] by including the following features: 

• The symbols should not be assumed to be independent of the space variables outside of a 
compact set. Otherwise, we would face a lot of difficulties with the smoothing procedure in 
the Nash-Moser iteration. 

• The remainders in the calculus of [12] should be smoothing operators when they were merely 
bounded operators on L 2 with a small (0(7 _1 ), 7 large) norm in [12]. Moreover, we desire 
more systematic and easily applicable criteria than in [12] for determining the mapping 
properties of remainders. 

The techniques in [12] heavily use the fact that all symbols are either Fourier multipliers or they 
are independent of the space variables outside of a compact set. One major goal here is to get rid 
of this assumption. Furthermore our symbols do not satisfy the standard isotropic decay in the 
frequency variables (basically there is one direction in frequency space in which there is no decay). 
We thus adopt a different strategy that is based on the Caldcron-Vaillancourt Theorem and more 



In fact, the loss in [6] is a loss of one singular derivative dx +/3 dg/e which implies a very bad control with respect 

to e. 
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specifically Hwang's proof of this Theorem [9]. Our motivation for doing so is that we shall only rely 
on L°° bounds for pscudodiffcrcntial symbols, while the classical proofs seem inapplicable when the 
frequency decay fails. The situation is even worse here because some "expected" results on adjoints 
or products of singular pseudodiffcrcntial operators seem not to hold. For instance, asymptotic 
expansions of symbols do not hold beyond the first term, and even the justification of the first term 
in the expansion depends on the order of the operators. Our final results are thus in some ways 
rather weak, but they seem to be more or less the best one can hope for in a singular calculus. 
Fortunately, the calculus is strong enough to be applicable to a variety of geometric optics problems 
for both wavetrains and pulses, including problems that display an amplification phenomenon. 

We thus review the results of [12] by improving them along the lines described above. For 
practical purposes, we have found it convenient to first prove general results on L 2 -boundedness 
of pscudodiffcrcntial and oscillatory integral operators. This will be done in Section 3 below. The 
calculus rules proved in Section 5 are then more or less "basic" applications of the general results. 
We have also found it convenient to include in the same article, the results for both the whole space 
and the periodic framework. Results in the case of the whole space are gathered in Sections 6, 7, 
8 and will be used in a future work to deal with pulse-like solutions to hyperbolic boundary value 
problems, while the companion article [G] is devoted to wavetrains. 
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Part A: singular pseudodifferential calculus for wavetrains 
2. Functional spaces 

In all this article, functions may be valued in C, <C N or even in the space of square matrices 
^#7v(C) (or C NxN ). Products have to be understood in the sense of matrices when the dimensions 
agree. If M G ^#jv(C), M* denotes the conjugate transpose of M. The norm of a vector x G C N 
is \x\ := (x* x) 1 ^ 2 . If x, y are two vectors in C N , we let x ■ y denote the quantity Y^j x j Vji which 
coincides with the usual scalar product in R N when x and y are real. 

2.1. Functional spaces on R d . The Schwartz space ^(R d ) of < t?°° functions with fast decay at 
infinity is equipped with the family of semi-norms: 

VJeN, Hulled), j := sup sup (l + \x\ 2 ) J / 2 \d«u(x)\. 

a6N d ,|a|<J i£l d 

When equipped with this topology, .y(R d ) is a Frechet space. We shall say that a sequence {uk)kez 
in ,y(R d ) has fast decay if for all polynomial P, the sequence (P(k) itfc)fcez is bounded in S*(M. d ). 
The Fourier transform on ,S fi {R d ) is defined by 

V/e J^(R d ), V£eR d , /(£):=/ e~ ix< f(x)dx. 



In particular, the Fourier transform is a continuous isomorphism on ,5^{R d ). It is extended to the 
space of temperate distributions J7"(W. d ) in the usual way. 

For s G R, we let H s (R d ) denote the Sobolev space 

H s (R d ) := {it G ,9>'{R d ) / (1 + \t\ 2 ) s/2 u G L 2 (R d ) j . 
It is equipped with the family of norms 

1 

7 :_ (2T) fl 

L 2 -norm on R d . 

write || • 1 1 o instead of || • ||o, 7 for the L 2 -norm on R d . For simplicity, we also write || • || s instead of 
,i for the standard iP-norm (when the parameter 7 equals 1). 



V 7 >1, VueH s (R d ), (Ml 2 , 7 :-tt— j / {-f + m'W^- 

{ Z7T > JR d 

The norm || • || 0j7 does not depend on 7 and coincides with the usual L 2 -norm on R d . We shall thus 



2.2. Functional spaces on R d x T. We now extend the previous definitions to functions that 
depend in a periodic way on an additional variable 9. We shall in some sense "interpolate" between 
Fourier transform and Fourier series. Let us begin with the definition of the Schwartz space. The 
Schwartz space ^(R^ x T) is the set of "if 00 functions / on R d x R, that are 1-periodic with respect 
to the last variable, and with fast decay at infinity in the first variable, that is: 

Va,^GN d , VjeN, f(x,e) G R d x R^x a d% <^/(x,0)) G L°°(R d x R) . 

Using the perdiocity of / with respect to its last argument 9, one can replace equivalently L°°(R d xR) 
by L°°(R d x [0, 1]). The Schwartz space y{R d x T) is equipped with the family of semi-norms 

VJGN, \\f\\^ xT) ,j := sup sup (1 + \x\ 2 ) J ' 2 \d a x 4f(x,6)\ . 

(a,j)GN d xN (x,e)eR d x[0S] 
\a\+j<J 

When equipped with this topology, y{R d x T) is a Frechet space. We let y'(R d x T) denote its 
topological dual, that is the set of continuous linear forms on y(R d x T). 

The "Fourier transform" on y(R d x T) is defined by considering Fourier series in 9 and Fourier 
transform in x. More precisely, we introduce the fc-th Fourier coefficient: 

V/ey^xT), VfceZ, Va-GR d , c k (f)(x) := [ c ~ 2 ^ ke f( x ,6) d9 . 



For all integer k, the Fourier coefficient ct(f) belongs to the standard Schwartz space y(R d ). We 
can therefore define its Fourier transform Ck(f). In all what follows, the sequence (ck{f))k£Z is called 
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the Fourier transform of /. When we only consider Fourier series in 9, we use the notation c k to 
denote the fc-th Fourier coefficient. When we only consider Fourier transform with respect to the 
first variable i £ R d , we use the classical "hat" notation introduced previously. 

The reader can check that the Fourier transform (c k (f))kez of a function / e <y(M d x T) is a 
sequence in J?^(R d ) with fast decay. The inverse Fourier transform is defined through the formula: 



f(x,9)=^2^- 1 (c k (f))(x)e 



2iirkS 



where j^ -1 stands for the inverse Fourier transform in ^(R d ). To summarize, the Fourier transform 
is an isomorphism between ,5 fi (R d x T) and the sequences {g k ) k ^z in ^(W 1 ) with fast decay. 

Let us now extend the Fourier transform to the set of temperate distributions <9"(M. d x T). For 
u G y"(R d x T), the Fourier coefficients c k (u) G y"(R d ) are defined by the formula 

VfcGZ, V.ge.^), {c k {ul 9 )^, m ^ m :=(„,«,(*) e- 2 '" * °) ^ (Rd xT) ^ {Rd xT) ■ 

It is straightforward to check that there exists a constant C and an integer J such that for all k G Z, 
there holds the continuity estimate 

y 9 e^(R d ), \(c k (u),g)^ m ^ m \ < C(l + \k\) J \\ 

The Fourier transform of u is the sequence {c k {u))k^i in ,5^"{R d ). For an appropriate constant that 
is still denoted C and a possibly larger integer that is still denoted J, there holds the continuity 
estimate 

VgE.y{W l ), |<^H,fl)^ (a a )i ^ (R a ) |<C(l + |fc|) J || fl ||^ (R1I ) )J . (1) 

The continuity estimate (f ) is uniform with respect to k e Z: the constant C and the integer J are 
independent of k. Moreover, the Fourier transform 

u G y'(R d x T) i — > (v&)) k & G ^'(R d f , 

is an isomorphism between 5^"{R d x T) and the sequences in ^"(W 1 ) that satisfy a uniform estimate 
of the type (1). The inverse Fourier transform is defined as follows: for a given sequence (itfc)fcgz 
in S"(R d ) satisfying a uniform continuity estimate with respect to k, we define an element v of 
y"(R d x T) by the formula 

V/ G ,y(M <i x T), (I'Jly^i^D.y^xT) : = ^2 (u-k^^ 1 ( c k(f))) ^, {Rd): y {Rd) , 

where J^" -1 denotes the inverse Fourier transform in S^(R d ). Indeed the reader can check first that 
v is well-defined, that it is a continuous linear form with respect to the topology of .^{W 1 x T) and 
that Ck{v) equals u k for all fceZ. 

For s G R, the Sobolev space H s (R d x T) is defined by 



H s (W l x T) := {u e y'(R d x T) / (c k (u)) k& 



G H s 



and 



(i +k 2 + \er\^)(o\ 2 ^<+^} 



V 7 >1, V U £^xT), || U || 2 >7 := * £ f tf + e + \tfy\c 



It is equipped with the family of norms 

- HT r + *~ + Kn" icfewuj' 2 

The norm |j • ||o )7 does not depend on 7 and coincides with the usual L 2 -norm on R d x T. We shall 
thus write || • || instead of || • || , 7 for the i 2 -norm on R d x T. More precisely, if / G L 2 (R d x T), 
then the Fourier coefficient 

c k (f)(x) := / e- 2i * k0 f(x,e)d6 
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is well-defined for almost every x G R d , and c k {f) belongs to L 2 (R d ) 
ity) . The Parseval-Bessel equality and Plancherel's Theorem give 



(use Cauchy-Schwarz inequal- 



\f(x,9)\ 2 dxd9 = Y, / \c k (f)(x)\ 2 dx=\\f\\ 2 . 

In what follows, we always identify the space L 2 (R d xT) and Fourier series in 9 G R whose coefficients 
belong to €2(^5 L 2 (R d )). For simplicity, we also write || ■ || s instead of || ■ || Sj i for the standard if s -norm 
on R d x T. 

Remark 1. Observe that our notation for the norm \\ ■ || Sj7 is consistent with the notation for 
functions that are defined on R d . More precisely, if u G H s (R d ), then one can also consider u as an 
element of H s (R d x T) that does not depend on 9, meaning that only the 0-th harmonic in 9 occurs 
(cq(u) — u and Ck{u) —Q ifk^O). The norms of u in H s (R d ) and H s (R d x T) coincide. This is 
the reason why we omit to write the underlying space R d or R d x T in the definition of the norms 



We now introduce the "singular" Sobolev spaces that we shall widely use in this article. From 
now on, we consider a vector f3 G R d \ {0} that is fixed once and for all. For set and e G ]0, 1], 
the anisotropic Sobolev space H s ' £ (R d x T) is defined by 



H s 



; (M d x T) := {u G y'(R d x T) /Vfc G Z, c k (u) G Lf oc (R d ) 



and 




2nk(3 



|clH(£)| 2 d£<+c»}. 



It is equipped with the family of norms 

V 7 >1, Vtie H s ' £ (R d x T) , 



2 



(2tt) c 



E 



1 2 + 



2irk/3 



\c k (um\ Z d^. 

(2) 

Let us observe that the definition of the space H s,e depends on e, and there is no obivous inclusion 
H s ' £l C H s ' £2 if £1 < £2 or £1 > £2- However, for a fixed £ > 0, the norms || • ||jj-«,e )7l and || ■ ||jja,s i72 
are equivalent. In particular, (2) defines a norm on the space H s > £ (R d x T) defined above. When 
to is an integer, the space H m ' £ (R d x T) coincides with the space of functions u g L 2 (R d x T) such 
that the derivatives, in the sense of distributions, 

01 / a \ a-d 



d Xl + — d g j ... (d Xd + -j- d g | 



belong to L 2 (R d x T) 
one derivative. 



In the definition of the norm 



ct\ + ■ ■ ■ + ad < to , 
e i7 , one power of 7 counts as much as 



In what follows, we shall also make use of the spaces ^ (I 



id 



x T), k G N: these are the spaces of 



continuous and bounded functions on t°xl that are 1-periodic with respect to their last argument, 
whose derivatives up to the order k exist, are continuous and bounded. 

3. The main L 2 continuity results 

Our goal is to develop in Section 4 a singular symbolic calculus on R d x T. This Section will give 
the basic results to achieve this goal. As in [12], the symbols that we shall consider do not satisfy 
the standard decay estimates in the frequency variable. Consequently, it will be more difficult to 
show that remainders in the symbolic calculus are smoothing operators. As a matter of fact, this 
property will hold only in a framework of the anisotropic Sobolev spaces defined above. A more 
cmbarassing consequence of this non-decay is that there seems to be little hope for developing a 
paradiffcrcntial version of the calculus below. More precisely, in the paradiffcrcntial calculus theory 
(see e.g. [11]), symbols have a fixed, say W k,ao , regularity in x. To cope with this small regularity, 
one introduces an isotropic frequency cut-off in the space variable. The regularized symbol belongs 
to the class S7\ and satisfies a suitable spectral condition, which yields continuity results for the 
associated pseudodifferential operator. This strategy applies when symbols are in £ with the 
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standard decay property (each derivative in £ yields one negative power of |£|). However, when this 
frequency decay does not hold or when it holds only in an anisoptropic way, the smoothing procedure 
yields symbols in the class Sq\ where derivatives in x are not balanced anylonger by derivatives in 
£. For such symbols, even with an appropriate spectral condition, there seems to be very little hope 
for continuity results in Sobolcv spaces. 

The remarks above are the main reason why we base our approach on the Caldcron-Vaillancourt 
Theorem [2] . More precisely, we shall prove continuity results on L 2 (R d x T) with symbols satisfying 
L°° bounds (no decay in the frequency variable is needed, see Theorem 1 below). This is the same 
strategy as in [12]. However we shall use more elaborate tools in order to get some refined estimates. 
Our goal is to get rid of the compact support assumptions in [12], and to lower the regularity required 
on the symbols whenever this is possible. We refer the reader to [4, 3, 1, 9] for some background on 
the Calderon-Vaillancourt Theorem and some generalizations. Here we clarify how these results can 
be adapted to a mixed situation where part of the space variables lie in M. d while the other space 
variables lie in the torus T. As far as we know, all previous versions were restricted to the case of 
R d or to the case of the torus. Our first continuity result is: 

Theorem 1. Let (7:M d xTxl <i xZ^ ^NxN ^ e a continuous function 2 that satisfies the property: 
for all a, (3 £ {0, l} d and for all j G {0, 1}, the derivative (in the sense of distributions) d™ d 3 g a 
belongs to L°°(E. d xTxI d xZ). 

For u G y(W l x T; C N ), let us define 

V(M)£M d xT, Op(<r)u(z,0) := V / e ix < e 2i7rfee a(x, 6, k) c^)^) d£ . 

Then Op(cr) extends as a continuous operator on L 2 (M. d x TjC^). More precisely, there exists a 
numerical constant C , that only depends on d and N, such that for all u G y{R d x T;C N ), there 
holds 

||Op(<r)u|| < C|H|| |M|o, with I a I := sup sup d^d j g d p ,a . (3) 



a^G{o,i} d je{o,i} 



L°°(K d xTxR d xZ) 



The proof of Theorem 1 below is analogous to the proof of [9, Theorem 2]. We emphasize that 
in the assumptions on the symbol a, no finite difference with respect to the index k G Z appears. 
This is in sharp contrast with for instance the paradiffcrential calculus on the torus developed in [8] . 
The fact that we do not need to estimate finite differences in k will be helpful in Section 4 when we 
consider singular pscudodiffcrcntial operators. 

Proof of Theorem 1. The proof of Theorem 1 combines two ingredients. First, the main estimate 
(3) holds when a is smooth, say e ^ ? °°, with compact support in all variables. Second, it is possible to 
approximate a symbol a satisfying the assumptions of Theorem 1 by a sequence (cr p ) p6 N of smooth 
symbols with sup p l\a p \\ controled by |||cr|||. The corresponding pseudodifferential operators Op(cr p ) 
converge in a weak sense towards Op(er). 

For smooth symbols with compact support, integration by parts and derivation under the integral 
show that Op(cr) u belongs to y(R d X T) if u does. In particular, Op(cr) u belongs to L 2 (R d x T). 
This integrability property is not so clear under the general assumptions of Theorem 1. 

Let us state more precisely our first point. 

Lemma 1. Let a G ^§°{W l x T x R d x Z; C NxN ), that is: 

(i) cr(-, ■, -, k) =0 except for a finite number of integers k, 

(ii) <r(-, •, •, k) is a < ^°° function on M. d x T x M. d for all k G Z, with compact support in its first 
and third variables. 

Then for all u G ^(R d x T; C N ), Op(cr) u belongs to y(R d x T; C N ) and the estimate (3) holds with 
a numerical constant C that is independent of a and u. 



Here and in all what follows, a function u that is defined on ff X Z is said to be continuous if for all k S Z, <r(-, k) 
is continuous on ff. The set will represent either Mr, or M. d X T or analogous sets. We adopt the same convention 
for differentiability properties. 
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Proof of Lemma 1. We make the proof in the case N = 1. When a takes its values in the space 
of matrices C NxN , the result applies for each component. Following [9], it is sufficient to prove an 
estimate of the form 



Op(cr) u(x, 9) v(x, 9) dxd9 



! xT 



<CMiM| |M|o, 



(4) 



for all u, v e y(Wl d x T; C). We define a function ip on R d by the formula: 

d 

\/y = ( yi ,... 7 y d )eR d , <p(y) := + ty^ 1 . 

3=1 

In particular, tp belongs to L 2 (M. d ). Applying Fubini's Theorem, we have 



Op(cr) u(x, 9) v(x, 9) dxd8 



! xT 



(2 7T) 



o~(x, 0, ^, fc) Ck(u)(y) v(x, 0) Ax d6 dy d£ . 



xT> 



Starting from the relation 



3 * (x-y)-£ 



and integrating by parts, we obtain 



3=1 



(2 7T) 



xTxI 



<rtt(x, 0, £, fc) t/(a;, £, A;) v(x, 0) dx d8 d£ . 



where we have used the notation 

d 

cjj : 



3=1 

We use the expression 



d p 

T(l-%)a, U(x,£,k):= e- i ve<p(x-y)c k (u)(y)dy. 



v{x,e) 



(2 



ce(v)(r)) dr] 



in (5) and apply Fubini's Theorem again. Then we use the relation 
and integrate by parts. These operations yield 



i= E E* 

ae{o,i} d ,je{o,i} q'<q 



E 



(5) 



e l ^e 2 "" a s fa, 5, 6 fc) <9"' £/ (ar, f , fc) ^(x, 9, f , jfc) dx d0 d£ , (6) 



! xTxI 



where the * coefficients denote some harmless numerical constants that only depend on a, a' ,j, and 
where we have used the notation 



(fc + 



Cf(u)(7?)d77 e 



,2i7r£0 



The result of Lemma 1 follows by applying Cauchy-Schwarz inequality to each integral in (6) (here 
the integral also includes the sum with respect to the index fc e Z). Each derivative <9" _Q dg a$ that 
appears in the right-hand side of (6) can be estimated in L°°-norm by a harmless constant times 
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the quantity \\a\\ defined in (3). We thus get (here and from now on, C denotes a positive numerical 
constant that may vary from line to line) 



\iV < c 



IHI 2 I E E / \d«U(x,t;,k)\ 2 dxdA V / \V(x,9,tk)\ 2 dxd6d^. 

Each term on the right-hand side is estimated by using the Parseval-Bessel equality and Plancherel's 
Theorem (see [9] for the case of M. d , here the incorporation of the additional periodic variable is 
almost straightforward) : 

£ / \dZU(x,t,k)\ 2 dxdt<CY,\\ck(u)\\ 2 <C\\u\\ 2 , 
E/ \V(x,9^,k)\ 2 dxd9d^<C\\v\\ 2 . 

fcgZ JK d xTxR d 

The proof of Lemma 1 is thus complete. □ 
To complete the proof of Theorem 1, it is sufficient to prove the following approximation result: 

Lemma 2. Let a : W. d x T x M. d x Z — > C NxN satisfy the assumptions of Theorem 1. Then there 
exists a sequence (a p ) pefi in ^ °°(R d x T x R d x Z:C NxN ) such that: 

(i) sup pgN 1 0p I < C I c I for some numerical constant C that does not depend on a, 

(ii) for all u,v € y(R d x T), there holds 



lim / Op(a p )u(x,6)v(x,0)dxd0 = / Op(cr) u{x, 0) v(x, 9) dx d6 . 

The proof of Lemma 2 follows by the classical truncation-convolution argument. We leave the 
details to the reader. The convergence property (ii) follows from the dominated convergence Theo- 
rem. 

Combining Lemma 1 and Lemma 2. we obtain the main estimate (4) not only for smooth symbols 
with compact support but also for the more general class of symbols that satisfy the assumptions of 
Theorem 1. In particular, the Riesz Theorem shows that Op(rr) u coincides almost-everywhere with 
an element of L 2 (M. d x T), and the conclusion of Theorem 1 follows. □ 

It is useful to observe that in the proof of Lemma 1 above, we do not need the symbol a to have 
compact support with respect to the space variable x. As a matter of fact, compact support with 
respect to the dual variables (£, k) is sufficient to justify all the calculations. This observation will 
be used in the proof of Lemma 3 below. 

Of course, the classical version of the Calderon-Vaillancourt Theorem in IR d now appears as a 
particular case of Theorem 1 (apply Theorem 1 with a symbol a containing only the 0-harmonic and 
that is independent of 6 and similar test functions u) , see [3, page 18] and [9] . In the proof of Theorem 
1, no finite difference with respect to k appears because there is no need to gain intcgrability for 
the function U with respect to the variable 9 (because the torus has finite measure). An even more 
direct explaination is that for a bounded sequence, the iterative finite differences are also bounded 
so the assumption would be redundant. 

Applying formally Fubini's Theorem to the definition of Op(a) u, we have 

O p (a)u(x,6) = -^ J2 [ e^ x -y^e 2 ^ k ^a(x,9^,k)u(y,co)dCdyd^. (7) 

The latter formula is rigorous e.g. for smooth symbols with compact support in (£, k). In order to 
prepare the results of symbolic calculus, our next goal is to obtain L 2 continuity results for oscillatory 
integral operators as in (7) with more general amplitudes a; namely we should allow a to depend on 
(x, 9) but also on the additional variables (y, w), see e.g. [12, page 144]. Our second main continuity 
result is the following: 
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Theorem 2. Let a :R d xT xR d xT xR d xZ ^ C NxN be a continuous function that satisfies the 
property: for all a, ft G {0, l} d , for all j, I £ {0, 1} and for all v G {0, 1, 2} d , the derivative (in the 
sense of distributions) d% d J s dl d% a belongs to L°°{R d x T x R d x T x R d x Z). Let X i G ^o°(K) 
and X 2 G ^ °°(R d ) sate/y Xl (0) = * a (0) = 1. 

XTien /or aZZ it £ J^(R d x T), the sequence of functions (Ts)s>o defined on R d x T by 



^ ' fceZ JK d xR d xT 

X 2 (5 £) cr(x, 0, y, w, £, k) u(y, w) d£ dy duj , (8) 

converges in ,y'(R d x T), as <5 tends to 0, towards a distribution Op(a)u G L 2 (R d x T). TTizs Zimii 
zs independent of the truncation functions xijX2- Moreover, there exists a numerical constant C, 
that only depends on d and N , such that there holds 

0p((7) U < C I CT I Amp || U || , 



with I a I Amp := sup sup sup 

a, / 36{0,l}< i j,;e{0,l} i/G{0,l,2} d 



■ (9) 

L°°(R ci xTxR d xTxR d xZ) 



The proof of Theorem 2 splits in several steps. The first point is to show that the conclusion 
holds for smooth symbols with compact support in (£, k). In this case, the convergence of the 
oscillatory integral as d tends to follows from the dominated convergence Theorem, and the proof 
of the continuity estimate (9) relies on some arguments that are similar to those used in the proof of 
Lemma 1. This first part of the proof of Theorem 2 is achieved in Lemma 3 below. The end of the 
proof of Theorem 2 consists in justifying the convergence of the oscillatory integral for an arbitrary 
amplitude and in showing that (9) still holds. This part of the proof relies on a regularization process 
as for Lemma 2. 

The process used in Theorem 2 that consists in introducing cut-off functions in the frequency 
variables and in passing to the limit will be systematically used in what follows in order to define 
oscillatory integral operators and to show some properties on such operators. To highlight the 
difference between standard pscudodiffcrcntial operators and oscillatory integral operators (for which 
the integrals do not converge in a classical sense), we always use the notation Op for oscillatory 
integral operators. In that case, the representation by a convergent integral only takes place when 
the amplitude is integrable with respect to the frequency variables (for instance, when it has compact 
support with respect to these variables). 

Proof of Theorem 2. We begin with the following generalization of Lemma 1. 

Lemma 3. Let a G c €^° (R d xTxR d xTxR d x Z;C NxN ) have compact support with respect to 
that is, there exists an integer Kq and a positive number Rq such that o~(x,6,y,u>,£,k) = 
as long as \k\ > Kq or |£| > Rq. 

Then all the conclusions of Theorem 2 hold and the oscillatory integral Op(er) u coincides with 
the function 

(x,6) G R d x — V / e l< - x -^<e 2lnk< - e - u K{x,e,y 7 uj,£„k)u{y,uj)d^dyduj. 

' ' fceZ "'R <i xR ci xT 

Proof of Lemma 3. Our strategy follows closely the proof of Lemma 1. In particular, we keep the 
same notation for the function ip on R d , and we make the proof in the case N = 1 for simplicity. 

First of all, since a has compact support in (£, k) and is bounded, the sequence (Ts)s>o defined 
by (8) is bounded in L°°{R d x T). Moreover, the dominated convergence Theorem shows that T$ 
converges pointwise, as 5 tends to 0, towards 

T(x, 6) := 7 -L 3 W e* C*-»K e 2 k a{x, 0,y,u, £, k) u{y,u) d£ dy dw . 

There is no ambiguity in the definition of the latter integral since the function to be integrated has 
compact support in £ and fast decay at infinity in y (the sum with respect to k only involves finitely 
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many terms). Applying again the dominated convergence Theorem, (Tg)s >0 converges towards T 
not only pointwise but also in J?"(M. x T). It thus only remains to estimate the function T in L 2 
in order to complete the proof of Lemma 3. We emphasize that the proof below does not assume 
compact support of a in x or y, which will be useful in Section 4. 
For v £ .Y(R d x T; C), let us define the integral 



T(x,0)v(x,O) dx d9 



! xT 



= TTU E / e^ x -y^e 2 ^ k ^a(x,e,y,^^k)u(y^)v(x,0)d^dxdedydio, 

{ Z7r > fcgZ JK d xE d xTxK' 1 xT 

where we have applied Fubini's Theorem. We first expand v as a Fourier series in 9: 

v(s,0) e 8 **", 
fez 

apply Fubini's Theorem, and integrate by parts with respect to 9 using the relation 

2 47 r(fc+£)0 _ 1 /-, , B s 2in(k+e)6 

We apply a similar manipulation for u, and we obtain the relation 

I = J2 I e * (X ~ VH e2i7Tk (e ~ W) ^(x,9,y,uj,^k) u(y, ui, k) v(x, 0, k) d£_ da; d6 dy duj , 

fe£Z JR d xR £i xTxR £i xT 

(10) 

where we have introduced the notation 
cr a := (1 - d e ) (1 -d^a, 

fez v ' fez v ' 

The latter manipulations are justified by the fact that both sequences (ce (it))^gz and (cf (^))fez have 
fast decay in j"(R d ). 

Let us now transform the expression of / in (10) by integrating by parts with respect to £. More 
precisely, we use the relation 

d 

e *(*-vK =tp{x-y) % + %) 2 e * , 
integrate by parts with respect to £ in (10) and obtain 



'xfxTxfxT 

x ^ (as, f , y, £, k) ip{x - y) u(y, u,k)ip(x- y) v(x, 9, k) d£ dx d9 dy duj , (12) 
where we have used the notation 

d d 

^ ■■= ^ = IK 1 - a - w - • 

3=1 3=1 

A crucial observation for what follows is that the new term <p(x — y) 2 in (12) yields integrability 
with respect to either x or y. 

Now we follow the argument already used in the proof of Lemma 1. We use Fourier's inversion 
formula, and write 

1 



v(x,e,k) = -— , / e l ^v(r,,e,k)drj, 

where for each k £ Z, the partial Fourier transform v(-, •, k) with respect to x belongs to the Schwartz 
space ,5^{R d x T). Then we apply Fubini's Theorem in (12), and integrate by parts with respect to 
x. As observed above, applying Fubini's Theorem has been made possible thanks to the new factor 



12 JEAN-FRANgOIS COULOMBEL, OLIVIER GUES & MARK WILLIAMS 

ip{x — y) 2 which makes the integral in x converge. We make the symmetric operation with u instead 
of v and integrate by parts with respect to y. Eventually, we obtain a formula of the form 



/= E E E *E 

o,/3G{0 : l} d a'+a"<a 0'+P"</3 fcGZ 



g i (x-J/)-£ e 2ink (0-u) 



l x& d xTxR d xT 

x tf-W a b (x, 9, y, u, £, fc) 

x {ti^+Ftpip - y)) U(y, u, f , fc) {d^' +fi "<p(x - y)) V(x, 9, f , fc) d£ da; d0 dy dw , (13) 

where the * coefficients only depend on <i, a, a', a", /?, /?', /?", and where we have used the notation 

1 

(27)° 



^ IS) jRd 

Let us now observe that each derivative d^~ a ~ a dy~^ ~^ &\> that appears in (13) can be bounded 
in L°°-norm by C ||| cr||| Amp , where the quantity ||| cr||| Amp is dchned in (9). Eventually, we apply the 
Cauchy-Schwarz inequality on (R d x T ) 2 x R d x Z in (13), and we thus need to estimate integrals of 
the form 

E / \d a ' + P'p(x-y)\ 2 \U(y,uj,Z 7 k)\ 2 d£,dxd6dyduj 7 

k< - z JM d x R d x T x R d x T 

and symmetric expressions in V . The latter integral is computed by first integrating with respect 
to (x, 9). Then we apply Plancherel's Theorem for transforming the integral in y into an integral in 
rj. Applying Fubini's Theorem, we can get rid of the integral in £ (see the above definition of U in 
terms of u) and we are left with estimating a quantity of the form 



E / |w(y>w,fc)| 2 d?/da;. 



feez 

where u is defined by (11). The latter quantity is estimated by using Parseval-Bessel's equality and 
Fubini's Theorem again. Eventually, we obtain 

E / \d a ' + ^^(x-y)\ 2 \U(y,iu,^k)\ 2 d^dxd9dyduj<C\\u\\l, 

k£Z JM d xTSL d xTxM d xT 

and a similar estimate holds for V in terms of ||u||o- We have thus proved that there exists a 
numerical constant C such that there holds 

|/| <C|H| A mpN|o|H|o- 

In particular, this yields the bound (9) when the amplitude a is smooth with compact support. □ 

Actually, the proof of Lemma 3 even shows the following stronger result which is encoded in the 
formula (13). 

Corollary 1. Let a E ^°°(M d x T x R d x T x R d x Z; C NxN ) have compact support with respect to 
(£, k). Let {Zi, . . . , Zm} denote the set of all derivatives dg d l u that appear in the definition 
(9) of the norm ||| ■ |||Amp- 

Then there exist some continuous bilinear mappings 

&L,...,Jf M ■ ^(K d x T) x y{R d x T) — ► L 1 (R d xTxl^xTxI^xZ), 
which are independent of a, that satisfy a continuity estimate of the form 

Vm = l,...,M, || = Sf m (u,-y)|| L i (K d )<TxR d xTxR ci xZ ) < C ||u|| \\v || , 
and such that for all u,v 6 ,5 fi {R d x T), there holds 

[ 6p(o-)u(x,6)v(v,6)dxd9 = V V ( (Z m o-)3? m (u,v)dxd9dydcjd£, (14) 

JR d xT TO=1 fcgZ JS. d xTxm d xTxM d 

where the expression of the function Op(cr) u is given in Lemma 3. 
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For a general amplitude a satisfying the assumptions of Theorem 2, we need to dchnc the limit, 
as 5 tends to 0, of the truncated oscillatory integrals (8). The goal is to show that formula (14), 
which holds for smooth amplitudes with compact support in (£, fc), also holds for the more general 
class of amplitudes satisfying the assumptions of Theorem 2. 

Let therefore a satisfy the assumptions of Theorem 2, and let us define the truncated amplitude 
as, S > 0, by 

a s (x, 6, y, u, £, fc) := xi(S k) x 2 (<5 £) a(x, 6, y, u, £, k) . 
The truncated amplitude as has as many derivatives as a in L°°. Moreover, there exists a constant 
C x that only depends on \ := (xi,X2) such that 

V5g]0,1], I a 5 1 Amp < C x I a ||| A m P • (15) 

Let us now consider a nonnegative function p G c tofi°(M. d ) with integral 1. We then define the 
regularizing kernels 

VnGN, p n (x) := (n + l) d p((n + 1) x) . 
We also consider the Fejer kernel 

,.(„:. ' ( "-«" + '»; 9 ' ) 2 , F„(0):=„ + 1 , 
n + 1 \ sm(7r 0) J 

that belongs to ^°°(T) and whose integral over T equals 1. Then we define the regularized amplitude 



(?8,n(%,Q,y,u,€,k) := / p n {x - x')F n (8 - 6')p n (y - y') F n (u - - £') 

jR d xTxR d xTxR t! 

cr^x', 0', y', w', fc) da;' d0' dy' dw' d£' . (16) 

It follows from the classical Theorems of calculus that for all 8 > and for all integer ngN, as, n 
belongs to ^ 6 °°(R d x T x R d x T x R d x Z; C WxAr ) and has compact support in (£, fc). Differentiating 
under the integral, we also have the bound 

Vn G N, I (T S ,n I Amp < IlkallAmp • (17) 
Moreover, since as is continuous, the sequence (c^rOneN converges pointwise towards as- 
For all u,v G ^(R d x T), let us define the integral 

7,5 := / T s (x, 6) v(x, 6) dxd6 , 

where the function Ts is defined by (8). Applying Fubini's Theorem, we have 

h = T^j V / e* e 2 4 - fc ^ (x, 0, y, u, f , fc) «(y, w) 0) dx dfl dy dw d£ . 

' fceZ ^R d xTxR d xTxR d 

(18) 

We also define the quantity Is. n that is the analogue of (18) with the amplitude as, n instead of as- 
The sequence (as, n )n&i is bounded in L 00 (R d x T x R d x T x R d x Z) and it is supported in a 
fixed compact set with respect to (£,fc). We can thus apply the dominated convergence Theorem 
and obtain that (I§ t n)neN converges towards Is as n tends to +oo. Moreover, we can apply Lemma 
3 to the amplitude a$ „ and derive the bound 



7,5, J 



xT 



Op(cr<5 jr i) u(x, 9) v(x, 6) dx d6 



< C |||a 5 ,„|||Amp ||«||o IMIo < C x lalllAmp |M| \\v\\ . 



where we have used (17) and (15). Passing to the limit as n tends to +oo, we obtain the uniform 
bound 

|7 6 | < C x |||(T||| Am p||u||o||w||o. (19) 
If we can prove that (Is)s>o has a limit as 5 tends to 0, and that this limit is independent of the 
truncation function \, then we shall have shown that the sequence of functions (Ts)s>o converges in 
y'(R d xT) towards some limit Op(cr) u. Moreover, the estimate (19) will show that this distribution 
coincides with a function in L 2 (M. d x T) satisfying (9). (If the limit of (Is)s>o is independent of x, 
then the constant in (9) is given by passing to the limit in (19) with one particular choice of x ) It 
therefore only remains to prove that (Is)s>o has a limit and that this limit is independent of x- 
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Since the amplitude ag n is smooth with compact support in (£, k), we can apply Corollary 1. Wc 
obtain that Ig n can be written under the form 

u r 

V/ {Z rn as, n )£' m {u,v)dxdedydLodZ. (20) 

m=l fca *xTx^xTxl^ 

We wish to pass to the limit in (20). We first observe that the derivative Z m as, n is obtained by 
differentiating under the integral sign in (16), that is 

Z m a s , n (x,6,y,u},£, k) = / p n (x - x') F„(0 - 6') p n (y - y') F n (ui - w')p n (£ - £') 

jR d xTxR d xTxR d 

Z m a s (x',6',y',u',g,k)dx' d0' dy' do/ d£' . 
Consequently, the right-hand side of (20) is a finite sum of terms that all have the form 

(q„ * h){v) f(v) dv, T := R d x T x R d x T x R d , 



with h G L°°(Y), / G £ 1 (T), and g n is the corresponding regularizing kernel. (Recall that the sum 
with respect to k in (20) involves finitely many terms, where the number of terms only depends on S 
and not on n.) Applying Fubini's Theorem, we can make the convolution kernel g n act on / rather 
than on h. This only replaces g n by g n with 

Qn(x,0,y,u,Q ■= Qn{~x, -0, -y, -u), -£) . 

Then we use the convergence of g n * f towards / in L 1 (this is a classical result of convolution that 
is unfortunately false in L°° and this is the reason why we need to switch the regularization kernel 
from one function to the other). Hence we can pass to the limit as n tends to +oo in (20), and 
obtain 

M . 

lim I Sn =I s =\2\2 (Z m a s )^ m (u,v)dxd6dydujd^. (21) 

In other words, we have extended formula (14) to the truncated amplitude as- 

It is now straightforward to pass to the limit as S tends to 0. Indeed each derivative Z m as can 
be decomposed under the form 

M 

Z m a = xi(S k) X2(8 Z m & + e m ,m' (5) Xi (<* k) X2,m,m' (5 £,) Z m > a , 

m' — l 

where X2,m,m' G '^^(W 1 ) and e m ,m'{S) tends to as S tends to 0. Wc can therefore apply the 
dominated convergence Theorem in (21), and obtain the expression 



m r 

lim Is = 2_. I (Zm c) (u, v) dx d6 dy duj d£ , 

m _i ,,_™ JR d xTxR d xTxM ,i 



m— 1 k£l 

from which it is clear that the limit is independent of x- The proof of Theorem 2 is complete. □ 

The proof of Theorem 2 even shows that the formula (14) still holds under the more general 
assumptions of Theorem 2, and that it actually defines the function Op(er) u G L 2 (M. d x T) in a 
unique way: 

Corollary 2. Let a : R d x T x R d x T x R d x Z ^ C NxN be a continuous function satisfying 
the differentiability assumptions of Theorem 2. Let the bilinear operators Jz? m , m = 1, ...,M be 
defined in Corollary 1. Then for all u G ^(R d x T; C N ) and for all v G ^(M. d x T; C), the function 
Op{a)u G L 2 (R d x T) satisfies (14). 

Remark 2. Let us assume now that in Theorem 2, the truncation functions XiiX2 do not necessarily 
satisfy Xi(0) = X2(0) = 1. Then the corresponding sequence of functions (Ts)s>o converges in 
y'(R d x T) towards xi(0) xa(0) Op(cr) it. 
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Let us observe that for an amplitude a that only depends on (x,0,£,k) and not on (y, w), then 
the oscillatory integral Op(cr) u coincides with Op(cr) u. This can be checked directly by applying 
Fubini's Theorem and the dominated convergence Theorem. In that case, the convergence of the 
sequence (Ts)s>o in 5^"{R d x T) is much easier to obtain. Since we shall use this argument in what 
follows, we state the result in a precise way. 

Proposition 1. Let er:R d xTxR d xTxR d xZ^ ^NxN ^ e a con ti nuous function that satisfies the 
differentiability assumptions of Theorem 2 and that is independent of its third and fourth variables: 
a(x,9,y,uj,t;,k) = a$(x,9,!;,k). Then for all u £ S^{R d x T), Op(cr) u coincides with the function 
Op(<Ttj)u defined in Theorem 1. 

For simplicity, a function defined on R d x T x R d x T x R d x Z that is independent of its third 
and fourth variables is equally considered as a function on R d x T x IR^ x Z, that is we use from 
now on the same notation for a and aj in Proposition 1. We hope that this does not create any 
confusion. The following result is a more precise comparison between oscillatory integral operators 
and pseudodifferential operators. It contains Proposition 1 as a special trivial case. It is also the 
starting point for the pseudodifferential calculus developed in the following section. 



Proposition 2. Let a £ ^°°(R d x T x R d x T x R d x Z; C NxN ) be an amplitude, and let the symbol 
a £ c (f§° (R d x T x R d x Z; C NxN ) be defined by 



Then the operator Op(er) — Op(er) coincides with Op(r) ; where the amplitude r £ ^' & 00 ( R d x T x R d x 
T x R d x Z;C NxN ) is decomposed as 

r(x,0,y,u,£,k) = rx(x,6,y,w,£, k) + R 2 (x,6,y,w,£, k + 1) - R 2 (x,9,y,uj,£,k) , 



We observe that the amplitude R 2 does not depend on y but it depends on u, so it does not enter 
the framework of Proposition 1 . 

Proof of Proposition 2. Let us first assume that the amplitude a has also compact support in (£, k). 
In that case, the symbol a has compact support in (£, k), and we can apply Proposition 1 and 
Lemma 3: 



a(x, 9, £, k) := a(x, 9, x, 9, £, k) . 



with 





P (ct) u(x, 9) - Op(cr) u(x, 9) = P (ct - a) u(x, 9) 




c i (x-vH e 2 i tt k y-u) (-(^ 0) V} W) ^ fc ) _ ~( a . ) 0] Xj W) ^ ^ ^ w j dy duj d ^ 




e * («-v).£ e 2 i * k (e-u) (-^ ^ X) W) ^ fc ) _ -(^ (9) aj) 0j ^ a?)) w) dy dw d£ . 



Let us start with the first term on the right-hand side. Applying Taylor's formula, we get 





Op(ri)u(a:,0). 
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All manipulations are made possible by the compact support assumption with respect to (£, k) and 
the fact that u belongs to J^(R d x T). 

Let us now study the second term in the decomposition of Op(a — a) u. By standard results of 
calculus, the function R 2 defined in Proposition 2 is 1-periodic with respect to 9 and w, and is smooth 
(namely, ^j, 00 ) with respect to all its arguments. (The reason why we divide by 1 — e -2»ir(a;— e) an( j 
not by w — 9 in the definition of R 2 is to keep the periodicity with respect to both 9 and u. However, 
this is of little consequence, and R 2 basically counts as one w-derivative of the amplitude a.) Wc 
apply Abel's transformation and obtain 

_L_ W e l <*-»K e 2 <* k («-«) (a{x, 9, x, u, k) - a(x, 9, x, 9, £, k)) u(y, u) dy <ko d£ 

= 7AlE / e^ x -yy^e 2 ^ k ^-e^^ k+1 ^ e -^)R 2 (x,9,u J ,^k)u(y,u J )dydu J ^ 
= 7T~u E / ° 4 ^""^ ° 2 1 nk ^ ( R *( x > 9 > w ' & fc+l)--R2(», 0, w, £, w) dy d^ d£ . 

^ ' fceZ -'R <i xTxR' 1 

We have thus proved the result announced in Proposition 2 under the additional assumption that 
the amplitude a has compact support in (£, k). 

When the amplitude a does not necessarily have compact support in (£, k), we approximate a by 
a sequence <r<5, 5 > 0, as in Theorem 2. We leave as an exercise to the reader the verification that 
for the corresponding sequence of amplitudes (rs)se ]o,i]j there holds 

Vwey^xT), lim 6p(r 5 ) it = Op(r) u , 

where the limit is understood in the sense of S^'iW 1 x T) (use Remark 2). This completes the proof 
of Proposition 2. □ 

We have only proved Proposition 2 for very smooth amplitudes. In the following Section, we shall 
extend this decomposition to amplitudes with finite regularity by the standard smoothing procedure. 
At this stage, we feel free to shorten some of the arguments in the proof when they are completely 
similar to what we have already explained. 

4. Singular pseudodifferential calculus I. Definition of operators and action on 

SOBOLEV SPACES 

4.1. Singular symbols and singular pseudodifferential operators. Following [12], we now 
introduce the singular symbols and their associated operators. The classes of symbols arc defined 
by first considering the following sets. 

Definition 1. Let q > 1, and let & C M g he an open set that contains the origin. Let m £ BL Then 
we let S m (<ff) denote the class of all functions a : 6 x R d x [1, +oo[-> C NxN such that 

(i) for all 7 > 1, er(-, -,7) is on G x M d , 

(ii) for all compact subset K of 6 , for all a £ N q and for all v £ N d , there exists a constant 
Ca.u,K satisfying 

sup sup sup ( 7 2 + |£| 2 )- (m - M)/2 IS? ^ 7)| < C a , v , K ■ 

v£K f£K d 7>1 

Let us now define the singular symbols. 

Definition 2 (Singular symbols). Let m £ TSL, and let n £ N. Then we let S™ denote the set of 
families of functions (<J e ,7)ee]0,l],7>i that are constructed as follows: 

V (x, 9, £, k) £ R d x T x R d x Z , a En (x 7 9, £, k) = a (e V(x, 9), £ + 2 ^ - , 7^ , (22) 

where a £ S m (ff), V belongs to the space ^"QR x T) and where furthermore V takes its values in a 
convex compact subset K of 6 that contains the origin (for instance K can be a closed ball centered 
round the origin). 



SINGULAR PSEUDODIFFERENTIAL CALCULUS 



17 



In Definition 2, we ask the function V to take its values in a convex compact subset K of & so 
that for all e s]0, 1], the function eV takes its values in the same convex compact set K. This 
property is used in several places below to derive uniform L°° bounds with respect to the small 
parameter e. 

For simplicity, we shall not mention that S 1 ™ depends on the open set &. (It will be convenient 
from time to time to let denote various possible open sets.) With a slight abuse in the terminology, 
we shall refer to the elements of S™ as symbols rather than as families of symbols. We hope that 
this does not create any confusion. 

To each symbol a = (a £ . 7 ) ee ]o,i] l7 >i € S 1 ™ given by the formula (22), we associate a singular 
pseudodifferential operator 0p e ' 7 (a), with e G]0, 1] and 7 > 1, whose action on a function u e 
y(R d x T; C N ) is defined by 

Op^ (a) u(x,8) -j^-f V / e"<e 2 ^ ke <j(eV(x,e)^ + ^^,j) c~^T)(Od£. (23) 



(2 7r) d ^ Jtl 

Let us briefly note that for the Fourier multiplier o~(v, £,7) = i £1, the corresponding singular operator 
is d Xl + ((3i/e)dg. We now wish to describe the action of singular pseudodifferential operators on 
Sobolev spaces. As can be expected from this simple example, the natural framework is provided 
by the spaces H s ' e defined in Section 2. The following result is a direct consequence of Theorem 1. 

Proposition 3. Let n > d + 1, and let a £ 5™ with m < 0. Then Op e ' 7 (a) in (23) defines a 
bounded operator on L 2 (R d x T): there exists a constant C > 0, that only depends on a and V in 
the representation (22), such that for all e G ]0, 1] and for all 7 > 1, there holds 

v«e y(R d x T) , ||Op £ <» u|| < Ho . 

Let us observe that if we compare Proposition 3 with [12, Proposition 1.1], we obtain the same 
result with slightly less regularity on V , and above all without the compact support assumption on 
the function V. The constant C in Proposition 3 depends uniformly on the compact set in which 
V takes its values and on the norm of V in < 4 d+1 . Even when we do not state it so clearly, all 
constants in the results below will depend uniformly on a finite number of derivatives of the symbols 
(or amplitudes). 

Proof of Proposition 3. We wish to apply Theorem 1, so the only thing to check is that the symbol 
a e , 7 defined by (22) satisfies a bound of the form 

C<j,v 



Vee]0,l], V 7 >1, |||a £ , 7 ||| < 



j\ m \ 



For instance, the proof of the L°° bound follows from Definitions 1 and 2. Let us recall that for all 
e G]0, 1], £ V takes its values in a fixed convex compact subset K C 6 (because K has been assumed 
to contain the origin, see Definition 2), so we have 

d\ mil 

< C ,o,k I 7 2 + 



tr[eV(x,0),t+ ^,7 



2irkB 




C 

< -r- r- 

The L°° bounds for the derivatives of a Ei7 follow by applying the Faa di Bruno formula for the 
composition of functions. We omit the details. □ 

Remark 3. The result of Proposition 3 does not rely on the scaling of the substitution in the 
representation (22). More precisely, the same result would hold with the substitution V(x,6) instead 
of eV(x,9). The only important point in the proof is the fact that the function substituted in the 
v-variable takes its values in a compact subset of 6 that is independent of e, and that sufficiently 
many of its derivatives belong to L°° . This fact will be used several times in what follows. 

There is no great difficulty in extending Proposition 3 to symbols of positive degree. 

Proposition 4. Let n > d + 1, and let a 6 S 1 ™ with m > 0. Then Op e,7 (a) defines a bounded 
operator from H m > e (R d x T) to L 2 (R d x T) with a norm that is independent of £,7: there exists a 
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constant C > 0, that only depends on a and V in the representation (22), such that for all e G]0, 1] 
and for all 7 > 1, there holds 

Vw e y{R d x T) , ||Op e ' 7 (a) u|| < C ||u||ffm.. i7 . 

Proof of Proposition 4- It is sufficient to write the symbol 



2nk(3 



-m/2 



7 



2nk[3 



m/2 



to observe that the symbol 

(x,6,£,k)\ — > a E ,y(x,0,£,k) 7 2 + 



2irk(3 



1/2 



belongs to S„, and eventually to observe that the Fourier multiplier with symbol 

2\ m/2 



7 2 + 



is an isometry from H m,£ (R d x T) - equipped with the norm || ■ ||jj m , Ej7 - to L 2 (M. d x T). 



□ 



The result of Proposition 3 can be made more precise when the degree m of the symbol is negative. 
We shall not deal with the general case m < since in what follows, the case m = — 1 will be our 
main concern. Our result is the following. 

Proposition 5. Let n > d + 2, and let a € S~ . Then Op e ' 7 (a) defines a bounded operator from 
L 2 (R d x T) to H 1 ' £ (M. d x T) with a norm that is independent of s, 7: there exists a constant C > 0, 
that only depends on a and V in the representation (22), such that for all e € ]0, 1] and for all 7 > 1, 
there holds 

V«e ,^(R d x T) , ||Op e '» u\\ H1 ^ n < C IMIo . 

Let us observe that the regularizing effect of Proposition 5 requires an additional space derivative 
on the symbol compared to the L 2 bound of Propositions 3 and 4. This is the first occurence in 
this article of the general principle that "symbolic calculus (and not only boundedness of operators) 
requires spatial regularity" . Here, we study the action of the composition 

d Xj + ^de\ Op e -». 

Proof of Proposition 5. We first observe that Proposition 3 already gives the estimate 

V« G ^{R d x T) , ||Op e ' 7 (a) u|L < - ||u|| . 

7 

Using the definition (2) of the norm || • || ff i,e 7 , we see that it only remains to prove some bounds of 
the form 



Vj = l,...,d, Vuey^xT), 



d Xj + ^d e 



<C\\u\\ . 



(24) 



with a constant C that is independent of e, 7, u. We prove such a bound in the case j = 1 (this is 
only to simplify the notation). 

We can differentiate under the integral sign in the definition of 0p e,7 (a) u, see (23), obtaining 

(d Xl +Y de ) Op en (a)u(x,9) = (T 1 +T 2 +T 3 )(x,8), 
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where we use the notation 



i x-E Jl i 7r k 6 

e 



d v *[eV(x,6),t+^^-, 7 )-ed Xl V(x, 



i x-£ J2 i it k ( 

e 



d v a[eV(x,9),£ + ^^,j) -d 6 V(x,6) 



Cfc(«)(Odf, 



fees 

The terms Ti and T2 fall into the framework of Proposition 3. Indeed, the function 

0T>(u,£,7) : = ^io"(«^)7) . 

belongs to S°, since <r belongs to S _1 . Consequently, the term T\ reads Op £n (a\,)u where the 
singular symbol a\, belongs to n > d + 2. In the same spirit, the term T 2 reads Op e ' 7 (aj)'u 
where the singular symbol aj belongs to n — 1 > d + 1 (use the substitution (e V, e W) with 

W := d Xl V in the symbol d v a(v,£, 7) • iu). We can thus apply Proposition 3 to estimate 7\ and T 2 
in L 2 (R d x T). 

The remaining term T3 does not fall directly into the framework of Proposition 3 since there is 
an e missing in front of dg V, so we do not exactly have a singular pseudodifferential operator as 
defined in (23). However, we can still apply Theorem 1 (see Remark 3) to the symbol 

(x, 9, £, k) .— ► d v a (e V(x, 9), £ + ^M, 7 ) . dg V(x, 9) . 



Since V belongs to ^ h n (K d x T) with n > d + 2, the latter symbol is bounded and it has exactly as 
many derivatives in L°° as required in order to apply Theorem 1, and the L°° bounds on the symbol 
are independent of e £ ]0, 1] and 7 > 1. We can therefore apply Theorem 1 in order to estimate T3 in 
L 2 (R d x T). The estimates of Ti, T 2 and T 3 yield (24), so the proof of Proposition 5 is complete. □ 

Remark 4. It would be tempting to extrapolate from Propositions 3 and 5 that symbols in S1 m , 
m G N and n sufficiently large, define pseudodifferential operators that act from L? to H m,E . This 
is true indeed, but unfortunately the operator norm seems to blow up with e as soon as m is larger 
than 2 (as soon as m is larger than 2, one faces a derivative (dg/s) 2 and the factor e~ 2 is too large 
when acting on the function eV). We thus need to pay special attention and check carefully each 
result one by one in order to prove uniform bounds. 

The proof of Proposition 5 can be adapted without any difficulty to show that singular pseudo- 
differential operators with symbols of degree act boundedly on H 1 ' 6 and not only on L? . We feel 
free to omit the proof of this result which will be useful later on. 

Lemma 4. Let n > d + 2, and let a£S[J. Then 0p e,7 (a) acts boundedly on /T 1 ' e (R d x T) with a 
norm that is independent of £,7.' there exists a constant C > 0, that only depends on a and V in 
the representation (22), such that for all e £ ]0, 1] and for all 7 > 1, there holds 

Vne y(R d x T) , ||Ot/<» u|| ff i.. i7 < C \\u\\^.. n . 

4.2. Singular amplitudes and singular oscillatory integral operators. The result of Propo- 
sition 3 can be generalized to singular amplitudes by using Theorem 2 instead of Theorem 1. More 
precisely, let us first define the classes of singular amplitudes. 

Definition 3 (Singular amplitudes). Let m £ R, and let n £ N. Then we let denote the set of 
families of functions (as,7) e e]o,i],7>i that are constructed as follows: 

V (x, 9, y, oj, £, k) e R d x T x R d x T x R d x Z , 

a^(x,6,y,u } ,£,k):=a(eV(x,e),eW(y,u } ),£+^^,j) , (25) 
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where a £ S m (^i x &2), V and W belong to the space x T), and where furthermore V , resp. 

W , takes its values in a convex compact subset K\, resp. Ki, of G\, resp. 61, that contains the 
origin. 

To each amplitude a = (S El7 ) e g]o,i],-y>i G A™ given by the formula (25), we wish to associate a 
singular oscillatory integral operator Op 7 (a), that would be defined (formally at first) by 

Vee]0,l], V 7 >1, Op'\a) := 6p(o e>7 ) , 



and the oscillatory integral operator Op is introduced in Theorem 2. The problem is that, at this 
point of the analysis, the operator Op has only been defined for bounded amplitudes that have 

— - e,7 

sufficiently many derivatives in L°°, see Theorem 2. We can therefore only define Op (a) for 
nonpositive degrees m. The following result generalizes [12, Proposition 2.2]. The proof follows 
exactly that of Proposition 3 above, except that we use Theorem 2 instead of Theorem 1. 

Proposition 6. Let n > d+ 1. and let a £ A™ with m < 0. Then for all e £ ]0, 1] and for all 7 > 1, 

the amplitude a Ei7 satisfies the assumptions of Theorem 2. Moreover Op (a) defines a bounded 
operator on L 2 (M. d x T): there exists a constant C > 0, that only depends on a, V and W in the 
representation (25), such that for all e £ ]0, 1] and for all 7 > 1, there holds 

C 



Vm £ y{w x T) 



- — e,7 

Op (a) u 



< 



7 



I II I 



Nlo. 



The derivatives d 3 g dy d u d^a en are computed in the classical sense and all of them are continous 
bounded functions on R d x T x R d x T x R d x Z. These derivatives are obtained by applying the 
Faa di Bruno formula. 

Remark 3 still applies, meaning that the result of Proposition 6 would still hold if we had made 
the substitution (v,w) — > (V(x, 6), W(y, ui)) instead of (v,w) — > (eV(x,9) 7 eW(y 7 uj)). Here, the 
small parameter e is not crucial in order to derive the uniform L°° bound on the symbol. 

In the same way as we proved a regularization effect for singular pseudodifferential operators 
with symbols of negative order, we are going to prove a regularization effect for singular oscillatory 
integrals operators when the amplitude has negative order and is sufficiently smooth. 

Proposition 7. Let n > d + 2, and let a £ A~ x . Then the oscillatory integral operator Op 7 (a) is 
bounded from L 2 (M. d x T) to H 1,£ (M. d x T). More precisely, there exists a constant C > 0, that only 
depends on a, V and W in the representation (25), such that for all e s]0, 1] and for all 7 > 1, 
there holds 

Vu £ Y(R d x T) , 



Op 7 (a) u 



< C\\u\\ . 



Moreover, the derivatives of Op 
sign. 



£.7 



(a) u are computed by differentiating formally under the integral 



Proof of Proposition 7. In order to prove Proposition 7, we need to go back to the definition of the 
oscillatory integral operator Op in Theorem 2. Let n > d + 2, a £ A^ 1 , and let u £ ,y(R d x T). 
Let now X i € ^ °°(R) and \2 G < jf oo (R d ) satisfy xi(0) = % 2 (0) = 1. According to Theorem 2, we 
know that the function 0p £ ' 7 (a) u £ L 2 (M. d x T) is the limit in y"(R d x T), as S tends to 0, of the 
sequence of functions 



T s (x,9) 



(2 



& i (x-y)-£ e 2i7r fe (8-ui) 



xR d xT 



£ V(x,6),EW(y,u 1 ),Z+^^-, 1 



Moreover, Proposition 6 already gives the estimate 



Op (a) u 



<-Ho- 
7 



(26) 
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Each function T$ is bounded. Moreover, we have shown in the proof of Theorem 2 that the 
sequence (Ts)se]o,i] is bounded in L 2 (R d x T) and converges in .y"(R d x T). This boundedness 
property follows from the relation (21), Corollary 1 and (15). 

Let j £ {l,..,,d}. We are going to prove that there exists a constant C, that is independent of 
S, s, 7, u, such that 



0, 



de)T s 



<C\\u\\ . 



(27) 



Combining (27) with (26), we shall obtain the result of Proposition 7. Indeed, the uniform bound 
(27) is sufficient to show that the limit of T s in y"(R d x T) belongs to H 1 ' 5 ^ x T). (Here, we use the 
classical weak convergence argument and the uniqueness of the limit in the sense of distributions.) 
We thus focus on the derivation of the bound (27) for j = 1. 

Each function Ts has ^f 1 regularity, and can be differentiated under the integral sign by applying 
standard rules of calculus. We obtain 

ft 



d Xl 



de\Ts = Ti ,« + T : 



2.5 



T-. 



3,5 : 



where, similarly to the proof of Proposition 5, we use the notation 
1 



T 1>s (x,0) :-- 



{2ir) d 



2 7rfc/3 1 



R d xR d xT 



6 + 



sV(x,8),eW(y,w),£ > + 



2irkf3 



J (a 



d v a 



! xT 



T 3 , S (x,6) :-- 



(2 7T) 



eV(x,0),eW{y,u),t 



27rfc/3 



->7 



y)-€ e 2iTfe(e- 



• £ 9 X1 V(a;,e) 



7 w(y,w)d^dydw. 



""(y, w )dC dy dw, 



! xT 



£ V(i l J) 1 £%w) 1 f + 



27rfc^ 



7 -a e y(x,0) 



d^ dydw . 



The singular amplitude appearing in the first term T\_$ belongs to , so we can apply the same 
argument as in Proposition 6 to estimate this term in L 2 (M. d x T). In the same way, the singular 
amplitude appearing in the second term T2,a belongs to A~_ lt so we can still apply the argument of 
Proposition 6. The amplitude appearing in the third term T^ t s does not fall into a representation of 
the form (25) because there is an e missing in front of dgV. However, we can still apply Theorem 2 
because this amplitude has sufficiently many derivatives in L°° and the L°° bounds are independent 
of e,7 (same argument as in Remark 3). The result of Proposition 7 follows. In particular, we have 
justified that the derivative 

(3 



d Xl + — d g j Op (a) u 

is computed by differentiating formally under the integral sign (meaning that the singular amplitudes 
that appear after differentiation yield well-defined oscillatory integral operators). □ 

- — • £,7 ^ 

Extending the definition of Op (a) to the case m > does not seem so clear at first sight. The 
trick of Proposition 4 does not apply anymore, and we need another argument that we detail now. 
Due to the application that we have in mind (see the companion article [(!]), we restrict to the case 
of amplitudes of degree 1, meaning that the growth at infinity is 0(|£| + |fc|). We do not claim 
that our criterion in Lemma 5 below is sharp. As a matter of fact, there is some hope that refined 
methods may yield a similar result with less regularity on the amplitude, but this is not our main 
concern here. We simply note that using sufficiently many derivatives to integrate by parts enables 
us to justify the convergence of the truncation process without the compact support assumption of 
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Lemma 5. Let a G A\, n > 3(d+l). Leixi G ^ °°( K ) a7ld X2 G tf§°(R d ) satisfy xi(0) = X2(0) = 1. 
T/ien /or all u G J^IR' 4 x T), £/ie sequence of functions (Ts)s>o defined by (8) £/ie amplitude 

a £j7 converges in ^"{M. x T), and i/ie Hmii is independent of the truncation functions Xi>X2- 

As in the case a G A^, we let Op 7 (a) denote the oscillatory integral operator associated with 
a G At this stage, this operator maps 5? into 5?' . 

Proof of Lemma 5. As in the proof of Theorem 2, our goal is to show that for all u, v G -y(M. d x T), 
the integral Is defined by 



T s (x,0)v(x,0) dx d6, 



xT 



with Tg defined by (8) (just replace the general amplitude a in (8) by a £j7 ), has a limit as 8 tends to 
0, and that the limit is independent of the truncation functions XiiXz- Applying Fubini's Theorem, 
we have (let us ignore from now on the powers of 2 tt that do not play any role): 



fcGZ 



e* x < c 2i7!ke X 2 {$ v(x, 0) U{x, 0, f , k) dx d0 d£ . 



xTxI 



with 



U(x,6,£,k) :-- 



R d xT 



as,-y(x, S,y,w, £, k) u(y, u) dy dw . 



We claim that it is sufficient to prove an estimate of the form 

icr(x,e^,fe)i<c( e ,7,o,«) T ^ J[ire' 

and the convergence of 1$ will follow from the dominated convergence Theorem (the constants 
may depend in a very bad way on s but this is no concern for us since we are only interested in 
the convergence of the integral for every fixed value of e). The L°° bound for U is obtained by 
multiplying by the factor 

d 

(l-2tnk) 3 IJ(l-»6) 3 , 
i=i 

and by integrating by parts. Observing that a G A\ with n > 3 (d+ 1), we claim that the amplitude 
a £j7 satisfies the following bounds for each fixed value of the parameters e, 7: 



\d?,d, 



(x,o,y,u,t,k)\ <c(i + |e| 2 + fc 2 ) 1/2 , \p\+e<3(d + l), 



y •-'ui "E,7 



and we thus get 



(l-2ink) 3 J|(l -i^) 3 U(x,0,^k) 
which gives the result. 



<C(l + |£| 2 + fc 2 ) 



2\l/2 



□ 



In the following paragraph, we shall see how the oscillatory integral operator defined in Lemma 
5 for amplitudes in A\ acts on singular Sobolev spaces. 

4.3. Comparison between singular oscillatory integrals operators and singular pseudo- 
differential operators. Theorem 3 below extends the result of [f 2, Proposition 2.3] to our frame- 
work in the case of bounded symbols, and is the main ingredient in Section 5 to prove the symbolic 
calculus results. 

Theorem 3. Lei a G A°, n > 2 (d + 1), be given by (25), and let a G 5° be defined by 

V (x, 0, £, k) G R d x T x R d x Z , a en (x, 0, f , k) := a (e V{x, 0), e W(x, 0), £ + <- 

Then there exists a constant C > such that for all e G ]0, 1] and for all 7 > 1, there holds 

VtiG y(M. d x T) , 0p £ ' 7 (S)M-0p £ ' 7 (a)u < — llullo ■ 

7 



(28) 
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If n > 2 d + 3, f/ien /or another constant C , there holds 

Vue^xT), 6p £ ' 7 (5)u-Op £ ' 7 (a)u <CM|tt|| 0) (29) 

uniformly in e and 7. 

Proof of Theorem 3. The proof relies mainly on Proposition 2, which gives the expression of the dif- 
ference Op (a) u — Op e,7 (a) u. As a matter of fact, Proposition 2 holds for very smooth amplitudes 
but using the standard regularization procedure, the result of Proposition 2 can be extended to 
amplitudes for which the remainder r defined in Proposition 2 satisfies the assumptions of Theorem 
2. In what follows, we are going to verify that under the assumptions of Theorem 3, the remainder 
r can be estimated in the norm ||| • |||Amp and we shall feel free to apply Proposition 2 in this finite 
regularity framework. 

Let us recall that the remainder r can be split as r = n + r-i with r\ also defined in Proposition 
2 and 7*2 is a finite difference in k (the amplitude r2 does not depend on y). Here we consider the 
amplitude a £j7 . We are first going to estimate the amplitude n, and then Eventually, we shall 
prove the regularization estimate (29). 

• The amplitude ri reads 

ri = T^ / d w <Tj Uv(x,9),eW((l-t)x + ty,w),£+—^-,ij ■ e d yj W{{\ - t) x + t y, u) dt , 

( 3 °) 

with dj := a £ S" 1 . To prove that Op(ri) is bounded on I? , we wish to apply Theorem 2 and 
we thus try to control ||| r"i ||| Amp - For instance, the L°° norm of r\ is estimated by using the decay of 
o~j with respect to the frequency variables and we obtain 

C £ 

\rx(x,9,y,u,£,k)\ < — . 

7 

When estimating derivatives, the worst case occurs when the derivative with respect to uj, the d 
derivatives with respect to x and the d derivatives with respect to y all act on the term d Vj W{{\ — 
t) x + ty,u). This requires having a bound for the 2 d + 2 first derivatives of W in L°°. Derivatives 
with respect to £ are harmless since they only add decay with respect to the frequency variables. 
Under the assumption of Theorem 3, we thus get a bound of the form 

Ce 

Fl Amp < , 

7 

which is even better than what we aimed at in (28). 

• The estimate of the term ri is more delicate and requires some attention. We first use the trick 
that appears repeatedly in [12], namely we write 

a e , 7 =a(e V(x, 9), 0, £ + + a s (e V(x, 9), e W(y, u), £ + ^^,7) • e W(y, u) , 

where crjj still belongs to S°. The first term on the right-hand side does not contribute to the 

— e,7 

difference Op (a) u — 0p e ' 7 (a) it, see Proposition 1. We can therefore focus on the second term for 
which we have an extra e. To avoid introducing some new notation, we still use a £j7 to denote the 
second term on the right-hand side. Then we have r2 = i?(-, k + 1) — i?(-, k) with 

a £ ^(x, 9, x, uj, £, k) - a en (x, 8, x, 9, £, k) 



R(x,6,uj,Z,k):=^ l i_ e -2i*(«-«) 

77— d u a e ~(x,6,x,6,£, k) , if u> 

2 l 7T 

Considering k as a real variable (and not only an element of Z) , there holds 
\r%{x, 8, ui,£, k)\ < sup \dkR(x, 9, uj, £, k + re) 

k€[0,1] 

< C sup \d u dk a E ^(x, 9, x, u>, k + re) | . 
9,w,«e[o,i] 



if cj ^ 6», 
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The fc-derivative of a £j7 introduces a 1/e factor times a frequency derivative of the symbol <7j . The 
1 j e factor is compensated by the e factor in the term e W and the frequency derivative of the symbol 
yields a decay of the form 

2\ -V2 



2(fc + K)7T^ 



We thus obtain a bound 



|r 2 (x,#,w,£,fc)| < - 

7 



The estimate of the derivatives of r 2 follow the same strategy. Here there is no y-derivative to control, 
and the worst case occurs when we take d derivatives in x, one derivative in 9, and one derivative 
in io. This requires having d + 3 derivatives of the functions V, W in L°°. Since d + 3 < 2 (d+ 1), 
we thus derive a bound of the form 

C 
< — . 



F2 I Amp 



7 



Combining with our estimate of n and applying Theorem 2, we already get (28). 

Before going on and proving (29), we make an important remark. In our estimate of r 2 , we have 
taken into account the finite difference with respect to k in order to make a frequency derivative 
appear, to the price of a 1/e but gaining a 1/7. We could have also estimated each term of r 2 , 
meaning the terms i?(-, k + 1) and i?(-, k), separately. If we had adopted such strategy, we would 
not have gained a 1/7 but there would have been no trouble with the 1 je term. More precisely, the 
amplitude r 2 satisfies a bound of the form 



II r 2 III Amp 



< Ce. 

• Our goal is now to prove (29). Following Proposition 7, the derivative 

ft 



d xi 



£ 



dg Op(ri) u 



is computed by differentiating under the integral sign provided that the amplitude has sufficiently 
many derivatives in L°°, and similarly for r 2 . We show how to estimate such derivatives under the 
assumption n > 2d + 3. Let us start with the terms involving r\, which are actually easier. There 
holds 



d xi + -j- dej Op(n)u = Op + ~~~^ ) ri J u + °P(^in) m + Op f-^- 9 e r! 

We recall that the amplitude r\ is given by (30). To control d Xl r\ in the norm ||| • |||Amp, one just 
needs an extra space derivative than in the previous step. The same argument holds for dgri. 
Consequently, under the assumption n > 2 d + 3, we get 

a 1 ^ C 11 11 
7 



Op^n) u + Op ( — !- 9 e ri 
In order to estimate the amplitude 



2ifc7r/3i 



ri 



we use the decomposition (30), where we recall that the o^'s belong to S _1 . Compared to the 
previous step, this amounts to working with the symbols i£i<7j , which belong to S°, and we thus 
get uniform L°° bounds in 0(e). Eventually, we have obtained the bound 



<cii 

.7 



Mlo ■ 



The remaining task is to control the analogous expression with the amplitude r 2 instead of r\ . 
To control the terms that involve d Xl ri or (/3i/e) dgr2, we use the above remark. More precisely, we 
estimate each term with R(-, k+ 1) and i?(-, k) separately, keeping the e factor to cancel the singular 
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term /?i/e. This requires only one more derivative on the functions V, W since we take one more X\ 
or 9 derivative of the amplitude. The most tricky term corresponds to 



Hi 



2iknf3 1 



1*2 



For this final term, we use the decomposition 



2 i k tt ft 



r2= [i £ 



2i(k+ 1)tt6i\ . , . ( ^ 2iknB 1 \ , IN 



2i7rft 



i?(-,fc + l). 



The last term R(-,k + l)/e has already been estimated at the previous step and satisfies an O(l) 
bound in the norm ||| • |||Amp- What remains is a finite difference in k which corresponds to the symbol 
i £i a instead of a (and then making the substitution with the singular frequency £ + 2 k tt j3/e). We 
apply the same strategy as in the previous step to make a frequency derivative appear, to the price 
of a 1/e. Since the ^-derivatives of i£i a belong to S° and thus satisfy uniform L°° bounds, we end 
up with the estimate 



d Xl + ~^ d e 



<C\\u\\ , 



which completes the proof of (29). 



□ 



The following result extends Theorem 3 to the case of amplitudes with degree 1. In particular, it 
will clarify the action of singular oscillatory integral operators on Sobolev spaces. 

Theorem 4. Let a G A} n , n > 3d + 4, be given by (25), and let a <E S„ be defined by 

V (x, 9, £, k) e R d x T x R d x Z , o C)7 (x, 9, f , k) := a \ e V(x, 9), e W(x, <?),£ + 2 ^ - , 

Then the operator Op 7 (a) — Op 6 ' 7 (a) is bounded on L 2 , namely there exists a constant C > suc/i 
i/iai /or aZZ e G ]0, 1] and /or 7 > 1, there holds 



Vu e J^(R d x T), Op (2) u- Op 6 ' 7 (a) u <C||u|| . 



(31) 



In particular, Op (a) maps H 1 ' 6 into L 2 and there exists a constant C > such that for all e G ]0, 1] 
and for all 7 > 1, there holds 



V« G ^(M d x T) , Op ' (a) it 



<C||«| 



HM l7 



Proo/ 0/ Theorem 4. Let u G J^(R d x T). Then we know that Op ' (a) it is the limit in y (R d x T), 
as J tends to 0, of the sequence (T$), with (ignore from now on the powers of 2 7r): 

T S (x, 9) := V xi (S k) f e* C*-»K e 2 - fe («— ) X2 (5 £) a En (x, 9, y, u, £, k) u(y, u) d£ dy da; . 

fegZ JR d xR d xT 

Using the result of Proposition 2 (with a finite regularity, which can be justified by the standard 
rcgularization procedure), we decompose as usual 



Ts = Ti.s + Op(n,j) u + Op(r 2 ,s) u , 



with 



T s (x,6) :=y2xi(6k) f e l ^< e 2 ^ k ^ X2 (<5£) a £ , 7 (x, 9, x, 9, £, k) u(y, u) d£dy duj , 

feez J**x* d xT 

and ri j( 5, T2,s are as in Proposition 2 but are obtained by considering the truncated amplitude 
Xi{Sk) X 2 (<5£)a £j7 . 
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It is easy to show that the sequence Tig converges in 5^' (and even in a stronger sense) towards 
Op e,7 (a) it, because one can first integrate in (y, ui) and use the decay of the Fourier transform of u. 
We are now going to compute the limit as 5 tends to of 

Op(ri^) u + Op(r 2 ,g) u . 
Using the general formula of Proposition 2, we have 

n,s = Xi( s k) X2{5 £) ri + - )> ]xi($k) 8^X2(8 / 8 y .a en (x,6,{l -t)x + ty,u,£,k)dt. 



X>(<5fc)%X 2 (<5£) l\ 
3=1 J ° 

1 d f 1 

n := - / 9 Vj d^a e ^(x,9, (1 - t)x + ty,u,^,k)dt . 
1 .1=1 Jo 



Since d^ j 'd En is a bounded amplitude, we can apply Theorem 2 for the convergence of the term 
Op(xi(<5 k) X2(<5£) ri) u. More precisely, we have c^c^a e A°_ x , re — 1 >3d + 3, and we therefore 
know that the limit of this term is Op(ri) u. Moreover, the operator Op(ri) acts boundedly on L 2 , 
uniformly in £,7. 

We now deal with the remaining term in n j. It is sufficient to prove that the singular oscillatory 
integral associated with the amplitude 

Xi(<5 k) 8^x2(8 / d yj a en {x,e, (1-t) x + ty,u,£,k) dt, 
Jo 

has a limit in 5?' as S tends to 0. Since d Vj a belongs to n — 1 > 3 d + 3, we can apply Lemma 

5 to this term. Together with the extra 5 factor, we have shown that the limit of Op(ri j a) u in 3" 
is Op(ri) u, and that this term is controlled in L 2 uniformly with respect to e, 7. 

The analogous term with r2,s is dealt with in a similar way. The finite difference with respect to 
k plays the role of the £ derivative and we can prove uniform bounds of the amplitude in the norm 
II • II Amp by using the same arguments as in the proof of Theorem 3. We feel free to skip the details. 
The action of Op 7 (a) on H l s is obtained by combining (31) with the result of Proposition 4 for 
Op e ' 7 (a). □ 

5. Singular pseudodifferential calculus II. Adjoints and products 

5.1. Adjoints of singular pseudodifferential operators. Our results on adjoints are very easy 
consequences of all the preliminary results in Section 4. Let us start with the case of bounded 
symbols. 

Proposition 8. Let a £ S^, n > 2 (d + 1), and let a* denote the conjugate transpose of the symbol 
a. Then 0p E ' 7 (a) and 0p e ' 7 (a*) act boundedly on L 2 and there exists a constant C > such that 
for all e G ]0, 1] and for all 7 > 1, there holds 

Vu G y(R d x T) , ||0p £ ' 7 (a)* u ~ 0p e < 7 (a*) u|| < - ||u|| ■ 

If n > 2d + 3, then for another constant C, there holds 

Vu G ,y(R d x T) , ||0p e ' 7 (a)* u - 0p e ' 7 (a*) «|| ffI .« i7 < C \\u\\ , 
uniformly in e and 7. 

Proof of Proposition 8. As in [12, proposition 2.4], it is sufficient to observe that if a en is defined by 

- — - £,7 — 

(22), the adjoint operator 0p e,7 (a)* coincides with the singular oscillatory integral operator Op (b) 
associated with the amplitude 

~ / 2 7T k (3 

b s ,j(x,0,y,uj,€,k) := a e>1 (y,uj,€, k)* = a £%,u),(+ — - — ,7 

Then we apply Theorem 3 and the conclusion follows. □ 
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Proposition 8 can be extended to symbols of degree 1 up to an additional regularity in the space 
variables (this high regularity is mainly required to give a precise meaning to oscillatory integral 
operators) . 

Proposition 9. Let a G S^, n > 3d + 4, and let a* denote the conjugate transpose of the symbol 
a. Then Op E ' 7 (a) and 0p £ ' 7 (a*) map H l e into L 2 and there exists a family of operators R 5,1 that 
satisfies 

• there exists a constant C > such that for all e €]0,1] and for all 7 > 1, there holds 

V« G y(M. d x T) , ||i? e ' 7 ii|| < C\\u\\ , 

• the following duality property holds 

\/u,v e y(R d x T), (Op £ - 7 (a)u,w) L 2 - (u,Op^(a» L2 = (i? £ < 7 u,v) L 2 . 

In particular, the adjoint 0p £ ' 7 (a)* for the L 2 scalar product maps H l ' e into L 2 . 

Proof of Proposition 9. The proof is quite similar to that of Proposition 8. First of all, the regularity 
assumption n > 3 (d + 1) allows to pass to the limit in the standard truncation process and to show 
that the adjoint (with respect to the L 2 scalar product) of the operator Op £,7 (a) coincides with the 
oscillatory integral operator Op (6) associated with the amplitude 

b e ,y(x,0,y,u,€,k) := a £i7 (y,w,£,fc)* = a 2 ^ - , 7 

Then we apply Theorem 4 and the conclusion follows. □ 

Let us observe that we only have proved a symbolic calculus " at the first order" , meaning that we 
have not proved that the adjoint operator Op £ ' 7 (a)* admits an asymptotic expansion with more and 
more smoothing operators. Even in the case of "t^ 00 regularity for the substituted function V, it is 
not so clear that the second order expansion holds with a uniformly bounded remainder in the scale 
of spaces H h ' e . This bad behavior is more or less the same as in Remark 4 (consider for instance 
the case of differential operators of order 2). 

5.2. Products of singular pseudodifferential operators. We still follow [12] and begin with a 
special case of products. 

Proposition 10. Let a, b G S®, n > 2 (d + 1). Then there exists a constant C > such that for all 
e G]0, 1] and for all 7 > 1, there holds 

V« G y{M. d x T), ||0p E < 7 (a)0p e ' 7 (6)*w-0p 6 < 7 (a&*H| n < - ||u|| . 

7 

If n > 2 d + 3, then for another constant C , there holds 

Vug ^(U d x T) , ||0p^ 7 (a) 0p e < 7 (&)* u - 0p £ ' 7 (a&*) u|| ffl ,. )7 < C \\u\\ , 
uniformly in e and 7. 

Let a G S„ , b G S® or a £ S„, b G S^, n > 3 d + 4. Then there exists a constant C > such that 
for all e G ]0, 1] and for all 7 > 1, there holds 

V« G y(R d x T) , ||Op £ ' 7 (a) Op e ' 7 (6)* u - Op £ ' 7 (a6*) u|| < C \\u\\ . 

Proof of Proposition 10. In each of the three possible cases, the main point is to observe that the 
operator Op e ' 7 (a) Op 6 ' 7 (6)* coincides with the oscillatory integral operator Op (c) associated with 
the amplitude 

Ce,j{x,8,y,uj,£,k) := a e ^(x,9,^,k) b E „(y, u, f, k)* . 

The result is well-known for amplitudes with a sufficient decay with respect to the frequencies, and 
it holds in a more general framework provided that all oscillatory integrals can be defined (which 
is the case under the regularity assumptions stated in Proposition 10). The conclusion then follows 
from either Theorem 3 or 4. □ 
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A main improvement with respect to [12] is that we can now deal with all kinds of products by the 
classical ** argument. This improvement has been made possible because we have already shown a 
smoothing property for some remainders in the calculus (compare with [12, Propositions 2.6, 2.7]). 

Proposition 11. Let a,b G S^, n > 2 (d+ 1). Then there exists a constant C > such that for all 
e G ]0, 1] and for all 7 > 1, there holds 

Vu G J^(K d x T), ||0 P e ' 7 (a)Op e < 7 (&)M- Op e ' 7 (a6) u|| < - ||u|| . 

7 

J/ n > 2 d + 3, i/ien /or another constant C , there holds 

Vit G x T) , ||Op e <» Op £ ' 7 (6) tt - Op e < 7 (a&) uH^i..^ < C ||u|| , 

uniformly in e and 7. 

Let a E S„,b E S® or a G S®,b G n > 3d + 4. XTien i/iere exists a constant C > suc/i i/iai 
/or e G ]0, 1] and /or 7 > 1, £/iere holds 

VuG J^(M <i x T) , ||0p e ' 7 (a) Op e ' 7 (6) u - 0p e ' 7 (a 6) u||„ < C \\u% . 

Proof of Proposition 11. Let us deal for instance with the case a, b G n > 2 0!+ 3. Then we have 

Op 6 ' 7 (a) Op e ' 7 (6) = Op £ ' 7 (a)0 P e ' 7 (6**) = Op 5 ' 7 (a) (o P £ ' 7 (6*)* + i?!: 7 ) , 

where we have applied Theorem 3 to the symbol 6* and denoted FflJ the smoothing remainder 
(mapping L 2 into H l e ). Thanks to Lemma 4, we know that Op e,7 (a) acts continuously on H 1 ' 5 , 
uniformly with respect to £,7, so the product Op e ' 7 (a) FciJ can be rewritten as a remainder of the 
form R^fJ. The product 0p e ' 7 (a) Op e,7 (6*)* is dealt with by applying Proposition 10. We end up 
with 

Op e ' 7 (a) Op e - 7 (6) = 0p £ < 7 (a b) + R^J . 
The only other interesting case is a £ S^, b £ S®, n > 3 d + 4. Then we write again 

0p £ ' 7 (o) Op £ < 7 (&) = 0p e ' 7 (o) (0p £ < 7 (&*)* + i?!; 7 ) , 

and we observe that the product 0p e ' 7 (a) R e _lJ acts boundedly on L 2 , uniformly with respect to e, 7 
(use Theorem 4). The product Op £ ' 7 (a) Op e,7 (6*)* is dealt with by applying again Proposition 10. 
We leave all remaining cases to the interested reader. □ 

A surprising fact is that the ** argument also applies for products of operators with degree —1 and 
1. We feel free to skip the proof that is entirely similar to that of Proposition 11. 

Proposition 12. Let a G S~ l ,b G S^, n > 3c? + 4. Then 0p e ' 7 (a) Op e,7 (6) defines a bounded 
operator on H ,£ and there exists a constant C > such that for all e G]0, 1] and for all 7 > 1, 
there holds 

VuG y(M. d x T) , ||0p £ ' 7 (a) Op £ ' 7 (o) u - Op e ' 7 (a b) u\\ H i,.„ < C \\u\\ Q . 

The analogue of Proposition 12 seems unfortunately untrue when the product is taken the other 
way round, meaning when the operator of order +1 acts on the left. This can be seen for instance 
by choosing for the left operator the singular derivative d x 1 + (f3\/e)dg. We are then reduced to 
showing a bound in H 1,s for the terms T2, T3 appearing in the proof of Proposition 5. Such a bound 
is available for T2 but not for the last term T3. This fact gives rise to a special treatment of +1, —1 
products in the companion article [6]. 

5.3. Garding's inequality. The exact same arguments as in [12, page 155] apply to prove the 
Garding inequality without any compact support assumption on the symbols. We just need slightly 
more regularity on the symbols in order to apply Propositions 8 and 11 above. 

Theorem 5. Let a G S° satisfy Reo~(i>, £, 7) > Ck > for all v in a compact subset K of 6 '. Let 
now a G Sq, n > 2d + 2 be given by (22), where V is valued in a convex compact subset K. Then 
for all 5 > 0, there exists 70 which depends uniformly on V , the constant Ck and S, such that for 
all 7 > 70 and all u G y(R d x T), there holds 

Re (0p e ' 7 (a) tt; u) L 2 > (C - 6) \\u\\ 2 . 
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5.4. Extended singular pseudodifferential calculus. Following [12, page 153], we can extend 
all the above results on boundedness/adjoints/products to the larger class eS m of functions a : 
& x R d x M 2d x [1, +oo[^ C NxN such that 

(i) for all 7 > 1, cr(-, ■, ■, 7) is ^on^x^x M. 2d , 

(ii) for all compact subset K of 6 ', for all a £ N q and for all v S N 3d , there exists a constant 
C a M,K satisfying 

sup sup sup( 7 2 + |e| 2 )" (m - kl)/2 |^^ c( T(«,e,C,7)| <C a , v , K . 

V ^ K (€,C)eR d xK 2d ,|£|<|c| 7>l 

For such symbols, we use the substitution v — > e V(x, 9), £ —> £ + 2 it k (3/e, C — > (£, 2 77 fc /3/e), which 
gives rise to extended singular pseudodifferential operators of the form 

We can also define extended singular amplitudes and compare the extended oscillatory integral 
operators with the above extended pseudodifferential operator. All results in Sections 4 and 5 are 
proved in the same way for this extended class because we have always relied on general boundedness 
result such as Theorem 1 or Theorem 2 and these results can handle symbols or amplitudes in the 
extended class. 

The main interest of defining this extended class is to be able to consider pseudodifferential 
cut-offs of the form 

(3 D 9 



x {eV(x,6),D x ,^,i 

where x is supported in the region \(i <C |£ 2 1 (here Ci is the placeholder for £ and (2 is the placeholder 
for 2 ir k (3/e). Such cut-offs are useful to microlocalize near the specific frequency (3. We again refer 
to [6] for further applications of these techniques. 
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Part B: singular pseudodifferential calculus for pulses 



In this second part, we briefly explain why all the results of the first part also give continuity 
and symbolic calculus results in the case where the additional space variable and associated singular 
frequency lie in R. There are some slight technical differences in the proofs and we pay specific 
attention to those terms that have the worst behavior with respect to the singular parameter e. 
As far as the notation is concerned, we feel free to use the same notation to denote new classes of 
symbols, amplitudes and so on, in order to highlight the similarities between Part A and Part B. 
We hope that this will not create any confusion. 

The variable in R d+1 is denoted (x, 9), x G R d , 9 G R, and the associated frequency is denoted 
(£, k). In this new context, the singular Sobolev spaces are defined as follows. We still consider a 
vector [3 G R d \ {0}. Then for s G R and e G]0, 1], the anisotropic Sobolev space H s ' e (R d+1 ) is 
defined by 



; (R d+1 ) := {u G ,y'(R d+1 ) /me L 2 oc (R d+1 ) 

and 



1 



e 



\u(£, fc)| 2 d£dfc < +oo| . 



Here u denotes the Fourier transform of u on R rf+1 . The space H 3 ' e (M. d+1 ) is equipped with the 
family of norms 

2\ s 



V 7 >1, Vue H s ' £ (R d+1 ), \\u\\'fi a ,e 



1 



(27T) 



d+1 



\u(Z,k)\ d£dfc. 



When m is an integer, the space H m ' e (M. d+1 ) coincides with the space of functions u G L 2 (M. d+1 ) 
such that the derivatives, in the sense of distributions, 



o Xd + — °t 



u , ot\ + ■ ■ ■ + a c i < m . 



belong to L 2 (M. d+1 ). In the definition of the norm || • \\h™' e .7, one power of 7 counts as much as one 
derivative. 



6. The main L 2 continuity results 

We adapt the Caldcron-Vaillancourt Theorems to our framework where symbols will enjoy a 
suitable decay property in the additional space variable 9. The decay will allow us to avoid requiring 
a control of fc-derivatives to prove L 2 -boundcdness. 



••NxN 



be a continuous function that satisfies the property: 

), where the 



id for all j G {0, 1} , (9) d% % d£ a belongs to L c 



td+i 



Theorem 6. Let a : 

for all a,f3 G {0, l} d 

derivative is understood in the sense of distributions and (•) denotes the Japanese bracket. 
For u G y(R d+1 ; C N ), let us define 

V(x,9)eU d+1 , Op{a)u(x,0) := 1 / e i e > <r(x, 9, £, k) k) d^ dfc . 

(2n) a + L J Rd +i 

Then Op(a) extends as a continuous operator on L 2 (R d+1 ; C N ). More precisely, there exists a 
numerical constant C, that only depends on d and N, such that for all u G J^(R d+1 ; C ), there 
holds 



||Op(o-) u|| < C|H|| |H|o, with Hal 



sup sup 

a,/3G{0,l} d ]£{0,1} 



L oo( R d+l x]R< i+l) 



There is very little to change compared to the proof of Theorem 1. More precisely, we can 
reproduce the proof of Theorem 1 by making the following modifications: first of all, we consider 
d+1 space variables instead of d, and forget about the additional periodic variable. Then we perform 
all integration by parts except the one with respect to the last frequency variables (which would 
be aimed at gaining a decaying factor in 9). In the end, we add a (9) weight to the derivatives 
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of the symbol a and use the square integrable (6*} _1 weight on d™ U(x,£,k). The final estimates 
(Cauchy-Schwarz and so on) are unchanged. 

The main point to keep in mind is that, in contrast with the standard Caldcron-Vaillancourt 
Theorem which requires controlling one ^-derivative in L°° , we can use instead a decay property in 
9 for the symbol and still prove L 2 boundcdness of the associated pscudodiffcrcntial operator. As 
we have already seen in the case of wavetrains, avoiding a control of ^-derivatives is crucial if we 
wish to prove some bounds that are uniform with respect to the singular parameter e. 

We now turn to the case of oscillatory integral operators associated with amplitudes defined on 
R d+1 x R d+1 x R d+1 . We feel free to skip the proof of the following Theorem which is again a slight 
adaptation (with similar modifications as described just above) of Theorem 2. 

Theorem 7. Let a : R d+1 x R d+1 x R d+1 ->■ C NxN be a continuous function that satisfies 
the property: for all a, (3 £ {0, l} d+1 , for all j £ {0,1} and for all v £ {0,1, 2} d , there holds 
( w ) d x.e d y,u d l dfctr belongs to L°°(R d+1 x R d+1 x R d+1 ). Let X 6 ^ °°(R d+1 ) satisfy x(0) = 1. 
Then for all u £ ^(R^ 1 ), the sequence of functions (Ts)s>o defined on R d+1 by 

T S (x, 9) := — i— T f e l «•(*-»)+* («—)) X (S £, S k) a(x, 9, y, u , £, k) u(y, u) d£ dfc dy dcj , 

(ZTT) a ^ 1 J R d+l xR d+l 

converges in 5^" (R d+1 ) , as 5 tends to 0, towards a distribution Op(er) u £ L 2 (R d+1 ). This limit is 
independent of the truncation function \- Moreover, there exists a numerical constant C , that only 
depends on d and N , such that there holds 



Op(o)u < C I o- 1 Amp IM|o , 



with I a I Amp := sup sup sup 

a,/3e{0,l} d+1 J'G{0,1} v£{0,l,2} d 



- 1 ) 



In Theorem 7, one exchanges a /c-derivative for a decay in to. Using a decay in 9 is also possible 
but it does not seem to be useful for our purpose. However, it will be useful below to have a version 
of Theorem 7 that combines both decay in u> and 9 and thus gets rid of all fc-derivatives. 

Theorem 8. Let a : R d+1 x R d+1 x R d+1 ->• C NxN be a continuous function that satisfies the 
property: for all a, (3 6 {0, l} d+1 and for all v 6 {0, 1, 2} d , there holds (9) (w) d% g d^ u d\ a belongs 
to L cc (R d+1 x R d+1 x R d+1 ). Let X & ^ oo (R d+1 ) satisfy x(0) = 1. 

Then for all u £ .^(R^ 1 ), the sequence of functions (Tg)s>o defined on R d+1 by 

T S (x, 9) := [ e* «•(*-»)+* («-«-)) x (5 £, S k) a(x, 9, y, w> £, k) u(y, lo) d£ dk dy dcu , 

converges in ^''(R d+1 ) 7 as 5 tends to 0, towards a distribution Op(er) u £ L 2 (R d+1 ). This limit is 
independent of the truncation function X- Moreover, there exists a numerical constant C , that only 
depends on d and N , such that there holds 

Op(er) u ^ < C I cr| Amp IM|o , 

with IHUmp := sup sup (u) dZ e d^ u d%a\\ . 

Of course, when the amplitude a in Theorems 7 and 8 does not depend on its third and fourth 
variables, we are reduced to the case of pseudodifferential operators. 

7. Singular pseudodifferential calculus I. Definition of operators and action on 

SOBOLEV SPACES 

Let us first define the singular symbols. 

Definition 4 (Singular symbols). Let m £ R, and let n £ N. Then we let 5™ denote the set of 
families of functions (a£, 7 ) e g]o,i],7>i that are constructed as follows: 

V (x, 9, £, k) £ R d+1 x R d+1 , a e „{x, 8, £, jfe) = a (s V(x, 6>), £ + — , 7 ) , (32) 
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where a £ S m (ff), (9) V belongs to the space c ^(M. d+1 ) and where furthermore V takes its values 
in a convex compact subset K of & that contains the origin (for instance K can be a closed ball 
centered round the origin). 

To each symbol a = (a e ,-y) e £]o,i],7>i £ given by the formula (32), we associate a singular 
pseudodifferential operator Op e ' 7 (a), with e e]0, 1] and 7 > 1, whose action on a function u £ 
^(R d+1 ;C JV ) is defined by 

Op=.-r( ) tt (M) : =75-4rT / ** i$ ' x+k0) ° U V{x, 9), £ + ^, 7) k) d£ dfc . (33) 

Let us briefly note that for the Fourier multiplier o~(v, £, 7) = the corresponding singular oper- 
ator is d Xl + (/3i/e) 9g. The main difference with respect to (23) is that now the singular frequency 
k/e lies in all R and not only in a discrete set. This modification will not be negligcablc in some 
places, especially when we compare the difference between oscillatory integral operators and pseudo- 
differential operators. Following Part A, we wish to describe the action of singular pseudodifferential 
operators on Sobolev spaces. The following result is a rather direct consequence of Theorem 6. 

Proposition 13. Let n > d + 1, and let a £ S™ with m < 0. Then 0p e ' 7 (a) in (33) defines a 
bounded operator on L 2 (R d+1 ); there exists a constant C > 0, that only depends on a and V in the 
representation (32), such that for all £ S]0, 1] and for all 7 > 1, there holds 

Vuey(R d+1 ), ||O P £ >»n|| <-^H| . 

The constant C in Proposition f3 depends uniformly on the compact set in which V takes its 
values and on the norm of V in ^ +1 . 

Proof of Proposition 13. The situation is more or less the same as in Proposition 3. The only thing 
to observe is that if the symbol a were independent of V, then we could not use Theorem 6 because 
{9} a en would not be bounded (except in the trivial case a = 0). We thus use the same trick as in 
[12] and write 

a £ , 7 = <7fo,e+— ,-y) (eV{x,6)^+—A -eV(x,9). 



Boundedness on L 2 for the Fourier multiplier is trivial, and as far as the second part is concerned, 
we can apply Theorem 3 (we could even apply the classical Calderon-Vaillancourt Theorem because 
taking only one fc-derivative is harmless for the second term since it contains an e factor). We obtain 
a bound of the form 

Vee]0,l], V 7 >1, ||K 7 -a(o,£+^,7)|||<^^. 

□ 

Let us again observe that the scaling e V is not crucial here to obtain a uniform L 2 bound, as 
in the case of wavetrains. This is because we still do not need to control any fc-derivative. The 
analogue of Proposition 4 works in exactly the same way. 

Proposition 14. Let n > d + 1, and let a £ S™ with m > 0. Then 0p e ' 7 (a) in (33) defines a 
bounded operator from H m,e (R d+1 ) to L 2 (M. d+1 ): there exists a constant C > 0, that only depends 
on a and V in the representation (32), such that for all s £ ]0, 1] and for all 7 > 1, there holds 

Vw G y(R d+1 ) , ||Op e ' 7 (aH| Q < C\\u\\ H m.,„. 

There is also a smoothing effect in the case m < that is analogous to the one proved in Proposition 
5 for wavetrains. 

Proposition 15. Let n > d + 2, and let a £ S~ x . Then 0p e,7 (a) in (33) defines a bounded operator 
from L 2 (M. d+1 ) to H 1,£ (M. d+1 ): there exists a constant C > 0, that only depends on a and V in the 
representation (32), such that for all £ E]0, 1] and for all 7 > 1, there holds 

Vu£ ^(R d+1 ) , ||0 P e < 7 (a) u\\ H1 , Sn < C Ho . 
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If one compares with the proof of Proposition 5, the analogue of the term T 3 does not have an 
additional e factor but depends linearly on dgV . Hence it falls into the framework of Theorem 6, 
uniformly in e and 7 (we use the decay in 9 rather than controlling one fc-derivative) . 

We can extend the above results to singular amplitudes which are defined in the following way. 

Definition 5 (Singular amplitudes). Let m G R, and let n£N. Then we let A™ denote the set of 
families of functions (a£. 7 ) 6 e]o,i],7>i that are constructed as follows: 

V (x, 9, y, u, £, k) G R d+1 x R d+1 x R d+1 , 

a En (x,6,y,u>,$,k):=a(EV(x,6),eW(y,u;),Z+^,"A , (34) 

where a G S m (^i x ff 2 ), (9) V and (u) W belong to the space ^ n (]R d+1 ), and where furthermore V, 
resp. W , takes its values in a convex compact subset K\, resp. K%, of &\, resp. &i, that contains 
the origin. 

Our continuity results of Propositions 6 and 7 extend to the case of pulses. For a G A r 7 " , our goal 
is to define, whenever the formula makes sense, the singular oscillatory integral operator acting on 
functions u G ^(IR^ 1 ) as follows: 



Ov'\a)u{x,6):=— \^ f e HH*-v)+H0-*>)) a (eV(x,6),eW(y,u),Z + ^,7 

u{y , id) d^ dfc dy dui . (35) 

We have the following result for bounded amplitudes. (The integral in (35) has to be understood as 
the limit in y'(M. d+1 ) of a truncation process in (£,&).) 

— -£,7 

Proposition 16. Let n > d + 1, and let a G A™ with m < 0. Then Op (a) in (35) defines a 

bounded operator on L 2 (R d+1 ): there exists a constant C > 0, that only depends on cr, V and W in 

the representation (34), such that for all e G ]0, 1] and for all 7 > 1, there holds 

Vue^(R d+1 ), Ov\a)u <-^Ho. 



7!' 



Proof of Proposition 16. We use the decomposition 



■ V(x, 9), 0, C + — , 7 )+ n U V ( x > °^ £ W ^ w). C + — » 7 ) ' e W(y, u) . 



Boundcdness on L 2 for the first term follows from Proposition 13 since we deal here with a singular 
pseudodifferential operator. As far as the second term is concerned, we can apply Theorem 7 because 
we have gained some decay in lo and taking one /c-derivative is harmless due to the e factor. □ 

The smoothing effect for amplitudes of degree —1 is more subtle. 

Proposition 17. Let n > d + 2, and let a G A^ 1 . Then Op ' (a) in (35) defines a bounded operator 
from L 2 (M d+1 ) into iJ 1,e (IR !i+1 ): there exists a constant C > 0, that only depends on a, V and W 
in the representation (34), such that for all e G]0, 1] and for all 7 > 1, there holds 

VuGJ^(M d+1 ), 6p e ' 7 (S)u <c|M| . 

Proof of Proposition 17. Proposition 16 already gives a control of 7 times the L 2 norm so it only 
remains to estimate the derivatives. Applying similar arguments as in the proof of Proposition 7, 

— - £,7 

the derivative (d Xl + Pi/ede) Op (a) u is computed by differentiating formally under the integral 
sign as long as the amplitude obtained by such differentiation defines a bounded operator onL 2 . We 
compute 

d xi + y de) Op £ '\a) u = T 1 +T 2 + T 3 , 
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where we use the notation 
1 



Ti(x, 



T 2 {x,0) :-- 



Tz{x, 



(2 






i 


( 




1 




(2 


lt) d + l 






d, 




Pi 





J (i-(x-y)+k(B~Lj)) 



i + 1 xK d + 1 



i (^i + ^T 1 ) ^£V(M),£%wU+7,7) u(y,u ) )d£dkdydu,, 



ft/3 
e 



(2 7T) 



d+1 



a « (C-(x-y)+fe(e-w)) 



ft/3 



9, a ( e 0), e W{y, uj)^ + ^, 1 )-d e V(x, i 



u(y,uj) d£dfcdy dw , 



t(y, w) d£ dfc dy dw . 



The L 2 bound of the terms Ti, T 2 follows from Propositions 16 since we deal with bounded amplitudes 
of sufficient smoothness. The only term to consider with care is T3. We decompose its amplitude as 
follows: 

k/3 



d v a[e V(x, 0), e W{y, w), £ + — , 7 • d e V(x, 9)=d v a[e V(x, 0), 0, £ 



kf3 



+ <r b [eV(x,8),eW(y,w), £+— 



n)-deV(x,6) 



■W(y,u>),d e V(x,e)] , 



where cr^f, iu, £,7) acts bilinearly on its arguments. The first term in the decomposition is a pseu- 
dodiffercntial symbol that has decay in and for which we can use Theorem 6 (the decay in implies 
that we do not need to control fc-derivatives). The second term in the decomposition has decay in 
both and u> and it even has an e factor. We can thus apply either Theorem 7 or 8 and derive an 
L 2 bound that is uniform in the parameters e, 7. □ 

The argument of Lemma 5 based on integration by parts still works for amplitudes of degree 1 , 
and we have 

Lemma 6. Let a 6 A^, n > 3 (d + 1). Then Op 7 (a) in (35) is well-defined from S fi (M d+1 ) into 
.J?"(W. d+1 ) as the limit of the operators associated with the amplitude x(<5£, S k) a, \ S ^0°° 

Singular oscillatory integral operators and singular pseudodifferential operators are closely linked. 
The results below are the direct extensions of Theorems 3 and 4. There is a new technical difficulty 
which arises because the set of ^-frequencies is no longer discrete and we thus really need to take 
derivatives while we had to deal with finite differences in Part A. However, the general ideas in the 
proof arc quite similar. 

Theorem 9. Let a G A° n , n > 2 (d + 1), be given by (34), and let a G 5° be defined by 

V (x, 0, £, k) G K d+1 x R d+1 , a en (x, 0, f , k) := a V(x, 9), e W(x, 0), £ + M 7 ^ . 
Then there exists a constant C > such that for all e G ]0, 1] and for all 7 > 1, there holds 
Vti £ ,y(M. d+1 ) , 6p £ ' 7 (5)u-Op e ' 7 (a)u 



< - No- 

7 



If n > 3d + 3, then for another constant C, there holds 



ViieyfO, Op"(a)u-Op en (a)u 



£,7, 



< c 



u\\o, 



uniformly in e and 7. 



The reason why we need 3 d + 3 derivatives on the symbol for the smoothing effect (rather than 
2 d + 3 as in Theorem 3) will be explained in the proof below. 
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Proof of Theorem 9. Let us first observe that when the amplitude a docs not depend on (y, w), there 
is no error in the difference Op (a) u — 0p e,7 (a) u, so we can restrict to the case where a has the 
form 

a E , 7 (x,6,y,LJ,t,k) := a leV{x,6),eW(y,w),£ + —,7 J ■e%w) , 

where cr(w, 7) • u>2 acts linearly on W2- 

• Following the ideas of Proposition 2, we can decompose the difference Op (a) — 0p e,7 (a) as 
Op {r\) + Op (r 2 ), with r x as in Proposition 2, and 

1 f 1 

r-2 ■= - / d u dka en [x, 6, x, (1 — t) 6 + tu>,£, k) dt . 
1 Jo 

The amplitude r\ reads 

n = i YlJ Q a i {eV{x,9),eW{{l-t)x + ty,uj)^+ ] ^-,^j ■ sd yj W((l - t) x + t y,oj) dt 

+ 7 J2j Aw(J 3 (eV(x,0),eW((l - t) x + ty,uj),^ + ■[sd yj W,sW}((l-t) x + ty,uj)dt , 

with o~j := a G S . To prove that Op^x) is bounded on L 2 , we wish to apply Theorem 7 and 
we thus try to control |||Amp- Since we can use some decay in w, it is sufficient to control one 
derivative in k, 2 d derivatives in £, d derivatives in x, d derivatives in y, one derivative in 9 and one 
derivatives in to. The worst case occurs when d derivatives in x and d derivatives in y act on the 
term d Vj W(. . . ) and we thus need W to have 2 d + 2 derivatives in L°° with the weight (lj). The 
factor e allows for a uniform control of the fc-derivative, and we thus get a bound of the form 

INlAmp <-, 

7 

with the quantity |||ti ||| Amp defined in Theorem 7. 

— ' £ '7 

• Let us now look at the operator Op (^2), which is more complicated. We compute 
7 ^ I &i a 3 (eV{x,9),eW{x,{l-t)9 + tuj),i+^- ll 



- V [ Pj a 3 (e V(x, 9), e W(x, (l-t)9 
7 E / & d *^ ( eV ( x > 



k B 

eW{x,{l-t)6 + tu),£ + —,'))-{ed bJ W 1 W](x,(l-t)6+tuj)dt, 



.1 

2-1 



with again (jj := a G S . The "leading" part of the amplitude r 2 is 

+ ^,7) -d u W{x,{\-t)e + tu)6t 



1 d r 1 ( 

~ E / Pj<rj(eV(x,0),0,t- 

1 o ( ,r, m n ,- k P \ W(X,6)-W(X,UJ) 

t E & °i [ s V ^ *)> o, e + 7J ■ f J e _ u 



for 9 ^ id. Using separate estimates for \6 — u\ < 1 or |0 — u\ > 1, we obtain 

C 

a,/3e{0,l} d+1 i>G{0,l,2} d v ; 7 

uniformly in the parameters e and 7. (Observe that we cannot take any fc-dcrivative because there 
is no e factor in the amplitude 7-3, but we have decay in two directions of M. 2 that are either (9 — lj, 9) 
or (9 — We are not exactly in the framework of Theorem 8 (where the decay takes place 

in the (0,uj) directions) but we claim that the continuity result of Theorem 8 still holds if one 
replaces the weight (9) (tu) by the above weight. The only important point is to have a weight in two 
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independent directions of the (9, a;)-plane. We can therefore conclude that the oscillatory integral 
operator Op(r3) is bounded on L 2 with an operator norm controlled by 1/7. 

• It remains to prove a bound in L 2 for the operator Op(r2 — r$), and r 2 — r% has the form 



1 d 



kj3 \ ( , fc/3 



ft la-j (eV(x,fl),eW(x > (l-t)e + tu;),C + -^,7j -o-j- I e ^(x, 0), 0, £ + -f, 7 



■9a,W(x,(l-t)0 + ta;)(i* 



3=1 

+ \Y1J & Aw(J J v ( x > e )> e W ( x > (1 - *) + * w), e + 7J ■ [e d u W, W] (x, (l-t)6+tu)dt. 
The L°° norm of this quantity is uniformly controlled by 

-el (6 + t (lj - 6)}-' 2 dt < -eiuj-6)- 1 . 
7 Jo 7 

The I/7 factor comes from the fact that aj belongs to S _1 , and the exponent —2 in the integrand 
comes from the fact that we have two functions that both have decay in their "fast" variable. Since 
we have an e factor available in the amplitude r 2 — 7-3, we can control one /c-derivative in L°° just 
using the weight (u — 6). Observe that we really need to control one fc-dcrivativc because this term 
has decay in one single direction of the ((9, w)-plane so we are not able to use the same argument as 
for the amplitude 7-3. 

At this stage, we have seen that each piece in the decomposition of the oscillatory integral operator 
Op (a) — 0p e ' 7 (a) gives rise to a bounded operator on L 2 with operator norm 0(1/7). We thus 
obtain the first part of Theorem 9. 

• In order to prove the smoothing effect, we need to control the first order singular derivatives 
( if the difference Op (a) u — Op e,7 (a) u. If we stick to the above decomposition n + r%, the piece 
r3 will be differentiated with respect to 9 and multiplied by 1/e. There will then be little chance 
to obtain a uniform control of this piece since we had no e factor there. We therefore use another 
decomposition of the amplitude and write 

- — . — Ej'y - — e - — £)T 

Op (a) u — 0p e,7 (a) u = Op (ri)u + Op (r2 ; j)u, 
where r\ is the same amplitude as above, and ri^ denotes the amplitude 

r 2 , t := <J (e V(x, 9), e W(x, u), Z + 7j-e W{x, w)-o (e V{x, 9), e W(x, 9), £ + 7^ -e W(x, 9) , 

(36) 

which each expression on the right-hand side has degree with respect to the frequencies (observe 
that here we have not applied Taylor's formula and integration by parts to get some decay in the 
frequency variables). 

The singular derivatives of the term Op (ri) u are computed according to the formula 

d Xl + 0p £ ' 7 (ri)w = Op 6 ' 7 U, Ki + ~~J u + Ov n {d Xl ri)u+^-Ov n {d g ri)u. 

Estimating each term in the above decomposition follows from arguments that were already used 
above. In particular, there is no problem here with the e factors since the 9 derivative on r-y yields 
an additional e factor, and it also yields some decay in the ^-direction. We can therefore prove 
a uniform L 2 bound for the singular derivatives of Op ' (r\ ) u as long as the regularity n of the 
functions V, W in the amplitude satisfies n > 2 d + 3 (compare with Theorem 3). 

e iT 

Let us now look at the singular derivative of the term Op (^2,0) u: 
(® Xl + 77") °P ' 7 ( r 2,t)) " = Op ' 7 f& + "~) r 2,j) u + Op' J (d Xl r 2 ,i)u+^Op J (d e r 2 .i)u. 

There is a subtletly here because the first amplitude on the right-hand side has degree +1 with 
repect to the frequencies, and this is the reason why we need n > 3 d + 3 in Theorem 9 (in order 
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to give a meaning to this quantity). For this first term, we use the Taylor formula and integrate by 
parts to get 



with 



J d u dk b E .j(x, 9, x, (1 — t) 9 + t u, £, k) dt^j u(y, uj) d£ dfc dy du; . 
[ ti + o [e V(x, 9), e W(y, w), £ + ^,7) • 



Since 6 Ei7 has degree +1, its fc-derivative has degree 0. More precisely, we can check that all terms 

arising when computing the derivative d u dk b £n yield an oscillatory integral operator that is bounded 

on L 2 for n > 3d + 3. 

— — 
The terms Op (d Xl T2,$) u and /3i/eOp (dgr2$)u are estimated by using the expression (36). 

Let us observe that the second term in the right-hand side of (36) is independent of (y,u>) so it 

gives rise to a genuine pscudodiffcrcntial operator (for which the continuity criterion of Theorem 

6 is less restrictive than the analogous result for oscillatory integral operators). Eventually, the 

interested reader can check that, using either Theorem 6 or Theorem 8, all amplitudes arising when 

computing d Xl r2,$ and fii/edgr^^ define oscillatory integral operators that are bounded on L 2 and 

whose operator norm is O(l) uniformly in £,7. We feel free at this stage to shorten the details that 

are very similar to many of the above arguments. □ 

In the same spirit as Theorem 4, we have the following result in the case of pulses. 

Theorem 10. Let a G A) i; n > 3 d + 4, be given by (34), and let a £ S„ be defined by 

V (ar, 0, C, k) 6 K d+1 x R d+1 , a S)7 (x, 0, £, k) := a [e V(x, 0),e W(x, 6), f + 7 
Then the operator Op (a) — 0p £ ' 7 (a) is bounded on L 2 , namely there exists a constant C > such 



that for all e € ]0, 1] and for all 7 > 1, there holds 

Vu G J^(R d+1 ), 6p^{a)u-0-p £n (a)u 



<CHk 



In particular, Op (a) maps H 1,£ into L 2 and there exists a constant C > such that for all e € ]0, 1] 
and for all 7 > 1, there holds 

Vtiey(l' f+1 ), 5p £ ' 7 (5W < CllttllHi..^. 



The proof is very similar to that of Theorem 4, with suitable modifications as in Theorem 9 in order 
to get some decay in the fast variables 9 and/or tu. 

8. Singular pseudodifferential calculus II. Adjoints and products 

The same results as in Section 5 hold in the context of pulses. We just state the corresponding 
results without proof in view of a future application to nonlinear geometric optics problems. The 
two first results deal with adjoints of singular pseudodifferential operators while the last two deal 
with products. 

Proposition 18. Let a £ 5^, n > 2 (d+1), and let a* denote the conjugate transpose of the symbol 
a. Then 0p £ ' 7 (a) and Op e ' 7 (a*) act boundedly on L 2 and there exists a constant C > such that 
for all e £ ]0, 1] and for all 7 > 1, there holds 

Vue.y(K. d+1 ), ||Op £ <»* U -Op £ - 7 (a*H| < - Ho. 

7 

If n > 3 d + 3, then for another constant C , there holds 

Vue ^(R d+1 ) , ||Op £ '»* u - Op £ ->*) u|| Hli . i7 < C ||u|| , 
uniformly in £ and 7. 
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Proposition 19. Let a G <5„, n > 3d+4, and Zet a* denote the conjugate transpose of the symbol 
a. Then Op 6 ' 7 (a) and Op e ' 7 (a*) map H 1,e into L 2 and there exists a family of operators R £r< that 
satisfies 

• there exists a constant C > such that for all e G ]0, 1] and for all 7 > 1, there holds 

Vuey(R d+1 ), \\R^u\\ <C\\u\\ , 

• the following duality property holds 

V«,ue y(R d+1 ) , (Op 6 ' 7 (a) u, v) L 2 - (u, Op e ' 7 (a*) v) L 2 = (i? e > 7 u, v) L 2 . 
In particular, the adjoint Op e ' 7 (a)* for the L 2 scalar product maps H 1 ^ into L 2 . 

Proposition 20. Let a,b G S^, n > 2 (d+ 1). Then there exists a constant C > such that for all 
£ £ ]0, 1] and for all 7 > 1, there holds 

VuG y{m. d+1 ) , ||0p^ 7 (a) 0p £ < 7 (&) u - Op £ < 7 (ao) u|| n < - ||u]] . 

7 

J/ n > 3 d + 3, i/ien /or another constant C , there holds 

Vu G y(R d+1 ) , ||0 P £ ' 7 (a) Op e ' 7 (6) u - 0p e ' 7 (a 6) u|| ffl ,. f7 < C \\u\\ Q , 
uniformly in e and 7. 

Let a G S„, b G S„ or a G S®, b € n > 3 d + 4. TTien i/iere exists a constant C > smc/i i/iai 
/or a// e G ]0, 1] and /or a// 7 > 1, there holds 

Vm G J^(R d+1 ) , ||0 P e - 7 (a) 0p £ ' 7 (&) u - Op £ ' 7 (a6) m|| < C ||u|| . 

Proposition 21. Let a G <S^\& G S„, n > 3d + 4. TTien 0p e ' 7 (a) Op e,7 (o) defines a bounded 
operator on Hl,s and there exists a constant C > such that for all e G]0, 1] and for all 7 > 1, 
there holds 

Vu G y(R d+1 ) , ||0p £ ' 7 (a) 0p e ' 7 (fe) u - 0p e ' 7 (a b) u\\ m , e ^ < C \\u\\ . 
Our final result is Garding's inequality. 

Theorem 11. Let a G S° satisfy Rcer(u,£,7) > Ck > /or all v in a compact subset K of (3 '. Let 
now a £ Sq , n > 2 d + 2 be given by (32), where V is valued in a convex compact subset K. Then 
for all 5 > 0, there exists 70 which depends uniformly on V , the constant Ck and S, such that for 
all 7 > 70 and all u G y(R d+1 ), th ere holds 

Re (Op E ^(a)u;u) L 2 > (C - S) \\u\\ 2 . 

There is of course an extended version of the singular calculus that allows for pseudodifferential 
cut-offs just as in the wavetrains case. 
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